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ABSTRACT 


The  presence  of  non-linear  behavior  in  potent iodynamic 
polarization  plots  has  resulted  in  difficulty  in  determining 
the  Tafel  constants  from  such  plots.  A  FORTRAN  based  program 
involving  numerical  differentiation  techniques  using  a 
graphical  display  Mas  used  to  determine  the  existence  of  the 
Tafel  regions. 

Various  alloys  polarized  in  synthetic  seawater  and  a 
3.5%  NaCl  solution  were  analyzed.  Although  severe 
concentration  polarization  often  dominated  the  cathodic 
branches  the  techniques  employed  did  allow  for  the  selection 
of  regions  which  approached  linear  behavior.  The  effects  of 
concentration  polarization  in  hindering  the  determination  of 
Tafel  constants  were  exemplified  by  the  uncovering  of  a 
cathodic  branch  containing  a  small  region  where  only 
actvivation  polarization  dominated  followed  by  the  onset  and 
total  domination  of  concentration  polarization. 

A  method  of  determining  where  the  anodic  and  cathodic 
currents  begin  to  dominate  the  potentiodynamic  polarization 


curve  is  introduced. 
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I-  INTRODUCTION 

Recent  studies  at  the  Naval  Postgraduate  School  have 
centered  on  the  vibrational  characteristics  of  certain  high- 
damping  alloys,  CPefs.  1-8].  Most  of  these  studies  have 
been  concerned  with  either  the  measurement  of  the  damping 
capacities  of  these  alloys  and/or  with  the  heat  treatments 
necessary  to  attain  the  microstructures  associated  with  the 
observed  high-damping  capacities.  In  addition,  two  separate 
studies  were  conducted  in  an  attempt  to  determine  the 
corrosion  characteristics  of  these  high  damping  alloys  in 
marine  environments. 

In  the  studies  involving  corrosion  characteristics,  both 
laboratory  and  actual  marine  environment  exposure 
experiments  were  conducted.  The  actual  marine  environment 
exposure  tests  calculated  corrosion  rates  based  on  direct 
weight  loss  methods.  The  laboratory  tests  conducted  were 
Linear  Polarization  and  Potentiodynamic  Polarization  in 
synthetic  seawater  [Ref.  9]  and  a  3.5%  NaCl  solution  [Ref. 
10]  environments.  Although  the  correlation  between  the 
laboratory  and  direct  exposure  methods  was  good,  the 
individuals  involved  in  both  of  these  laboratory  studies 
expressed  concern  with  the  difficulty  in  determining  the 
Tafel  constants  necessary  to  determine  corrosion  rates. 


The  actual  experiments  were  conducted  on  the  EG&G 
Princeton  Applied  Research  Model  351  corrosion  measurement 
system.  The  Modc^l  351  utilizes  a  PAR  Model  lOOO 
microcomputer  featuring  two  Motorola  68000  microprocessors 
and  a  touch-screen  input  to  create  an  environment  which  is 
well  suited  for  conducting  various  corrosion  tests.  The 
system  is,  however,  limited  by  the  lack  of  a  key  board  and 
more  importantly  by  an  operating  system  which  is  proprietary 
in  nature.  These  two  combine  to  make  any  analysis  of 
experimental  data,  other  than  by  the  means  provided  with  the 
operating  software,  impossible  unless  the  data  can  be 
transfered  to  another  computer. 


I I .  DATA  TRANSFER 

The  PAR  351  and  its  associated  softmare  has  the  ability 
both  to  plot  the  results  of  an  experiment  and  to  send  just 
the  data  points  to  a  printer  via  an  RS  232  serial  port. 
Although  the  software  supports  only  Houston  Instruments 
plotters,  the  printer  option  is  configurable  for  numerous 
printers  as  long  as  the  printer  supports  DCl  (XON,  ASCII  17) 
and  DC3  (XOFF,  ASCII  19)  protocol.  Other  RS  232 
requirements  such  as  baud  rate,  parity,  word  length,  and 
stop  bits  are  configurable  as  the  user  sees  fit. 

During  the  initial  data  transfers,  in  the  present  work, 
Hayes  Microcomputer  Inc.  SMARTCOM  II  modem  software  was 
used  to  capture  the  data.  This  method  was  successful  most 
of  the  time  but  since  the  two  systems  were  not  using  modems 
but  instead  were  directly  connected  via  their  serial  ports 
some  problems  were  experienced  and  editing  of  the  data  file 
after  transfer  was  always  necessary.  The  most  serious 
detraction  was  the  need  to  wade  through  various  menus  in 
order  to  accomplish  the  transfer. 

In  view  of  the  both  the  problems  experienced  with  using 
SMARTCOM  II  and  the  neccessity  for  the  user  to  be  familar 
with  the  way  SMARTCOM  II  works,  it  was  decided  to  write  a 


data  transfer/capture  program  which  would  both  speed 
execution  and  eliminate  the  need  for  user  knowledge  of 
communications  software. 

Using  the  sample  TTY  data  transfer  program  included  in 
the  manual  for  Microsoft  6W-BASIC  version  2  [Ref.  11]  as  a 
core  a  BASIC  program  was  written  to  accomplish  the  transfer. 
The  core  was  modified  to  eliminate  the  transfer  portion  and 
most  importantly  line  feeds  and  carriage  returns  sent  by  the 
3S1  were  automatically  deleted  and  a  carriage  return  was 
inserted  only  at  the  end  of  a  line  containing  data.  This 
eliminated  the  problem  of  blank  lines  appearing  within  the 


data  file  corresponding  to  page  headers  and  footers.  Other 
modifications  which  were  made  involved  checking  for  and 
correcting  two  formatting  variations  which  resulted  in 
compatibility  problems  between  the  transfered  file  and  its 
use  in  FORTRAN  application  programs.  The  first  involved  a 
change  from  a  floating  format  to  a  exponential  format  when 


the  potential  was  equal  to  0.0000  volts.  The  351  normally 
sends  potential  with  4  significant  digits.  When  the 
potential  is  0.0000  volts  it  is  displayed  as  -XXE-12.  The 


formatting  shift  is  not  compatible  with  any  FORTRAN  formats 
as  the  total  field  width  changes  in  addition  to  the  minimum 


number  of  available  digits.  This  occurance  is  now  trapped 
and  the  file  modified  to  reflect  the  actual  potential  of 
0.0000  volts.  The  second  formatting  variation  always 
occured  and  did  so  in  a  region  which  was  of  extreme  interest 
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as  follows.  When  the  current  density  changes  from  negative 
to  positive  there  is  a  range  where  the  exponent  is  E— 12. 


Unlike  FORTRAN  write  statements  which  always  reserve  a  space 
for  the  +  or  -  sign,  ie.  the  +  sign  is  not  printed  unless 
specified,  the  351  used  the  space  for  the  -  sign  if  present 


or  for  the  first  digit  if  not.  Again  this  resulted  in  a 
change  in  total  field  width  at  the  occurance  of  -a.bcdE-12. 
The  FORTRAN  read  statment  would  read  this  to  be  -a.bcdE-1. 
This  was  handled  by  simply  deleting  the  least  significant 
digit  such  that  the  number  is  read  as  -a.bcE-12  when  this 
occurs. 

In  its  final  form  the  data  transfer  program  requires  the 
user  to  reply  only  a  to  a  (R)eceive  or  (E)xit  prompt  and  to 
enter  a  file  name  if  the  (R)eceive  option  is  selected. 
Because  the  351  can  not  write  a  selected  configuration  to 
disk  the  user  is  still  required  to  properly  configure  the 
351  to  the  proper  settings  when  powering  up.  These  are: 

PROTOCOL  -  PRINT 

PARITY  -  NONE 

BAUD  RATE  -  2400 

STOP  BIT  -  1 

WORD  LENGTH  -  S 

Once  these  options  are  properly  set  on  the  351  the 
desired  data  file  should  be  loaded  into  the  351  and 


formatted  into  the  linear  option.  At  this  point  the  user 


selects  the  (R)eceive  option  on  the  Z248,  presses  the  “PRINT 


DATA"  prompt  on  the  screen  oi  the  351  and  enters  a  filename 
on  the  Z248. 

The  final  version  of  the  data  transfer  program  was 
compiled  using  Microsoft's  BASIC  compiler  Version  1.0.  By 
eliminating  the  writting  of  the  data  to  the  screen  in 
addition  to  the  file  and  using  the  compiled  program  vice  the 
BASIC  interpreter  a  file  which  will  use  60K  of  disk  space 
can  now  be  transferred  in  under  S  minutes  as  opposed  to  the 
20  minutes  it  took  with  SMARTCOM  II.  A  listing  of  the 
program  is  included  in  Appendix  B. 


ill.  CORROSION  RATES 

A.  THEORETICAL  EQUATIONS 

In  general  the  corrosion  behavior  of  a  metal  can  be 
described  by  the  use  of  the  Nernst  equation,  activation 
potential ,  concentration  potential  and  mixed  potential 
theory. 

The  Nernst  equation  is  a  modification  of  the  Gibb's 
Free  Energy  equation  and  expresses  the  tendency  of  a 
corrosion  reaction  to  proceed  in  terms  of  EMF  based  on  the 
activities  of  the  products  and  reactants  of  the  1/2  cell. 


E  =  E*:  -  ^log,(Q) 
nF 


E  -=  Electrode  oxidiation  potential 
(0  for  reduction) 

E'^  =  Equilibrium  oxidiation  potential 
(0  for  reduction) 

R  =  Gas  constant 
T  =  Temperature  ‘^K 

n  =  Number  of  electrons  taking  part 
F  =  The  Farday 

Q  =  Activities  of  products/activities  of 
reactants 

The  potential  of  the  cell  is  the  algebraic  sum  of 


oxidation  and  reduction  1/2  cells. 


Activation  polarization  describes  the  current  density 


which  exists  when  the  cell  is  polarized. 

-  ic.exp( +ZFaE/RT) 

=  i«,exp(-ZF(l~a)E/RT) 

io  =  equilibrium  exchange  current  density 
Z  =  Valance  of  the  active  species 
E  =  Polariztion  about  the  equilibrium  potential 
a  =  fraction  of  E  which  applies  to  i«nadi.c; 

R,T,F  as  above 

By  manipulating  the  above  equation  the  activation 
polarization  can  be  written  as: 

=  B  logxo  (i««/ie) 

0  =  2.3RT/aZF 

Although  activation  polarization  occurs  for  both  the 
anodic  and  cathodic  reactions  the  cathodic  reaction  can  also 
be  effected  by  concentration  polarization.  When  this  occurs 
the  reduction  reaction  occurs  at  such  a  rate  as  to  decrease 
the  concentration  of  '.ons  in  the  solution  near  the  surface 
of  the  metal.  This  results  in  a  change  in  the  potential  of 
the  1/2  cell  reaction  since  the  effective  activity  in  the 
Nernst  equation  changes. 
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The  concentration  polariztion  e  is  expressed  as: 


e  =  2.3RT  logi.0  (1  -ia^ii_) 
nF 


id  =  operating  current  density 
ii_  =  limiting  current  density 


The  limiting  current  density  is  inversely  proportional 
to  the  thickness  of  the  diffusion  layer.  This  thickness  is 
dependent  on  various  geometric  and  environmental  factors  and 
can  often  only  be  determined  experimentally.  The  effect  of 
concentration  polarization  can  be  minimized  by  agitating  the 
solution,  thus  decreasing  the  layer  thickness,  or,  as  in  the 
case  of  oxygen  reduction,  by  heating  the  solution  to  reduce 
the  concentration  of  the  dissolved  gas.  Both  of  these  means 
of  dealing  Mith  concentration  polarization  may  significantly 
alter  the  desired  corrosion  environment  and  distort  the 
correlation  betMeen  the  experimental  results  and  the  actual 
in-service  corrosion  rates. 

In  the  absence  of  concentration  polarization  a  dynamic 
polarization  diagram  of  a  cell  reaction  Mould  appear  as  in 
Figure  1.  The  slope  of  the  anodic  and  cathodic  curves 
d0/d(logxo  ic«)  would  be  given  by  B  as  described  in  the 
activation  polarization  expression. 

Most  cells  are  affected  by  concentration  polarization 
such  that  the  cathodic  polarization  '^  =  The 

effect  of  this  is  shown  in  Figure  2.  It  is  obvious  from 


s 

3* 
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this  that  lagj.0  (ia)  approaches  logio  (ii-)  the  linear 
behavior  of  the  polarization  curve  is  distorted  and  in  cases 
where  it_  and  io  are  not  sufficiently  separate,  a  linear 
region  may  not  exist. 

Mixed  potential  theory  concerns  concerns  the  immersion 
of  a  single  piece  of  metal  into  an  electrolyte  and  the 
subsequent  reactions.  The  theory  consists  of  two  simple 
hypotheses. 


1.  Any  electrochemical  ''raction  can  be  divided  into  two  or 
more  partial  oxidiation  and  reduction  reactions. 


2.  There  can  be  no  net  accumulation  of  electrical  charge 
during  an  electrochemical  reaction.  [Ref.  12  :  pp  314] 


It  is  these  hypotheses  that  allows  for  the 
determination  of  two  important  parameters  necessary  for 
estimating  corrosion  rates.  When  the  cathodic  branch  of  the 
1/2  cell  with  the  highest  equilibrium  potential  intersects 
the  anodic  branch  of  the  other  1/2  cell  the  corrosion 
potential  dcor-i-  and  the  corrosion  current  density  icar-r-  can 
be  used  to  estimate  corrosion  rates. 

When  the  corrosion  system  consists  of  only  a  single 
metal  in  contact  with  a  solution  which  contains  a  single 
reduction/oxidation  reaction  the  overall  reaction  is  the 
linear  sum  of  the  reduction  currents  and  the  oxidation 


currents . 


» 


If  the  solution  contains  more  than  one 
reduction/oxidation  reaction  the  overall  reaction  may 
proceed  as  the  linear  combination  provided  that  the 
equilibrium  exchange  current  densities  of  the  secondary 
reduction  reactions  are  within  one  order  of  magnitude  of  the 
main  reduction  reaction.  The  secondary  oxidiation  reactions 
proceed  in  the  same  manner.  As  a  result  the  secondary  redox 
reactions  affect  either,  both,  or  none  of  the  branches  of 
the  overall  reduc tion/oxidiation  reaction.  [Ref.  13:  pp. 

199] 

If  concentration  polarization  is  not  present  both  the 
anodic  and  cathodic  polarization  curves  would  consist  of  a 
of  horizontial  lines,  representing  the  current  added  by  the 
reaction  at  that  potential,  and  slanted  lines  of  some  slope 
0  representing  the  activation  polarization  of  all  the 
oxidation  or  reduction  reactions  which  are  influential  at 
that  potential.  In  the  presence  of  concentration 
polarization  the  possibility  of  introducing  numerous  non- 
linearities  into  the  cathodic  branch  is  distinct.  As  the 
overall  reduction  reaction  is  the  sum  of  all  the  influential 
reduction  currents  any  region  could  consist  of  zones  where 
the  activation  polarization  of  the  ith  reaction  dominates 
the  concentration  polarization  of  the  (i-l)th  reaction. 

This  can  make  the  potentiodynamic  curve  extremely  difficult 


f 
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B.  DETERMINATION  OF  CORROSION  RATES 

As  previously  described  the  presence  of  numerous 
electrodes  and/or  metals  in  the  corrosion  system  can 
significantly  alter  the  shape  of  the  polarization  curve. 
Despite  this  if  it  mas  possible  to  design  a  current 
measurement  system  that  mas  capable  of  measuring  only  the 
anodic  and  cathodic  currents  the  determination  of  icoi-i- 


could  still  be  simply  made  since  logi.o  ic 


Mould  be  at  the 


intersection  of  the  anodic  and  cathodic  branches  of  the 
polarization  diagram. 

In  the  real  Morld  polarization  curves  take  on  a  'Y' 


i 

i 


shape  as  a  result  measuring  instrument  limitations  and  the 
complexity  of  the  corrosion  environment.  can  be 

measured  directly  since  it  is  the  equilibrium  potential  of 
the  system.  As  the  potential  is  scanned  tomards  Pcor-.-  the 
net  current  (ianaaics  ~  ic:««rM»axe)  is  measured.  Only  Mhen 
the  potential  is  displaced  far  enough  from  does  one  of 


■s 

% 

»?■ 


the  tMo  currents  overwhelm  the  other  and  does  typical  Tafel 
behavior  take  place. 

The  presence  of  numerous  1/2  cell  reactions, 
concentration  polarization  and  the  need  for  one  of  the 


currents  to  overwhelm  the  other  can  all  combine  to  make  the 


determination  of  difficult. 

If  concentration  polarization  is  not  present  at  some  60 


from  0c;e>.'.-  the  polarization  diagram  should  adopt  true  linear 
behavior.  The  amount  of  60  at  which  the  curve  should  adopt 
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linear  behavior  is  in  itself  open  to  disscussion.  Ailor 
states  that  60  should  be  at  least  50mV  from  0ee>r-.-  [Ref.  14: 
pp.  199].  Once  the  linear  regions  are  found  two  methods  can 
be  employed  to  determine  ieor*-*  The  simplest  is  to  draw  two 
straight  lines  from  the  linear  regions  towards  0eof-t'>  The 
two  lines  should  intersect  at  ieoi-f  A  second  method  is  to 
measure  the  slopes  0  (Tafel  constants)  of  the  linear  regions 
of  the  anodic  and  cathodic  branches.  These  values  can  be 
used  in  conjunction  with  data  determined  from  a  polarization 
resistance  test  to  determine  the  value  of  ieor>i'  by  use  of 
the  equation: 

i-^ar-r-  — _ 0e0.« _ 

2.3(B.+B=)PR 

f 

If  concentration  piolarization  is  predominant  over  the 
cathodic  region  of  the  polarization  curve,  the  so-called 
'knee'  method  is  often  employed.  In  this  method  a  tangent  is 
draixn  on  the  knee  of  the  cathodic  branch  such  that  the 
angles  formed  on  both  sides  of  the  tangent  with  the  cathodic 
branch  are  equal.  When  this  has  been  accomplished  either  of 
the  methods  discussed  in  the  previous  paragraph  may  be  used. 

Many  other  methcsds  have  been  developed  to  determine 
icor-t-  or  to  determine  corrosion  rates  directly  from 
polarization  data.  Princeton  Applied  Research  uses  a  Chi- 
squared  minimization  technique.  This  technique  is  included 
in  the  software  package  which  operates  their  Model  351 
corrosion  measurement  system  as  PARCalc.  By  by  their  own 


The  logxo  (1  -  icj/ii_)  term  can  be  replaced  by  the  series 


expansion : 


loqio(l+X)  =  2.3C  X  -  X  =  ♦  X  =  -  X -•  .. 

2  3  4 

where  X  =  -i^/ic. 

The  final  equation  to  be  solved  would  be: 


0  -  0eo^r-  =  Bc:logio(i«a)-Belogio(ic»)+2.3[Z(-l)"-^*-.X"J 


The  difficulty  with  this  approach  is  that  a 
sufficient  number  of  terms  must  be  carried  in  the  series  in 
order  to  accurately  represent  the  logj.o  term.  This  causes 
two  severe  problems.  The  first  is  n->-2  data  points  must  be 
used  to  solve  the  system  of  equations.  When  one  considers 
that  i««  is  normally  on  the  order  of  uA/cm^,  exponentiating 
it  to  the  nth  degree  can  easily  result  in  a  matrix  to  be 
salved  involving  coefficients  of  SO  or  60  orders  of 
magnitude.  The  end  result  is  that  while  the  system  can  be 
salved  such  that  the  coefficients  determined  will  result  in 
the  accurate  representation  of  the  actual  overvoltage,  the 
coefficients  themselves  have  been  influenced  by  round-off 
error  and  catastropic  cancellation,  and  represent  merely 
numbers . 

The  other  problem  exists  with  attempting  a  solution  of 
this  form.  As  discussed  previously  with  the  possibility  of 
mixed  electrodes,  the  associated  redox  processes  may  or  may 
not  effect  the  overall  behavior  of  the  system.  Even  if  the 


limiting  current  densities  ii_  were  all  known,  it  would  still 
be  necessary  to  determine  their  effect  on  the  overall 
process  in  order  to  include  them  in  the  system  of  equations 
to  be  solved.  As  a  result,  for  m  different  electrodes,  m 
experiments  involving  Just  the  electrode  and  the  metal  would 
have  to  be  carried  out  Just  to  determine  their  influence. 
Consider  the  use  of  such  a  method  with  synthetic  seawater. 

At  least  10  different  chemical  substances  could  effect  the 
overall  reaction  in  addition  to  the  major  metal  oxidiation 
and  hydrogen  reduction  reactions.  It  should  be  obvious  that 
if  the  limiting  current  densities  are  not  known  the  size  of 
the  system  of  equations  to  be  solved  could  easily  exceed  the 
number  of  available  data  points.  More  importantly,  one 
would  have  to  assess  the  advantages  of  performing  a 
polarization  in  a  mixed  medium  if  the  same  test  has  to  be 
performed  separately  on  the  individual  components  first. 

C.  CALCULATION  OF  TAFEL  CONSTANTS 

In  order  to  determine  the  location,  if  one  existed,  of 
the  anodic  and  cathodic  Tafel  regions  two  numerical 
dif ferentation  techniques  were  employed.  These  two  methods 
were  the  Four  Point  Central  Difference  method  and  the  ur »  of 
a  Cubic  Spline  interpolating  polynomial.  A  graphical 
display  of  the  derivatives  allowed  the  user  to  determine  the 


v.-vj  vi  v:  wj  v^.-  v.-  v.-yjwy 


w 


region  of  linearity.  Once  this  region  was  determined  a 
first  order  linear  regression  is  performed  on  the  original 
data  within  the  region  to  determine  the  Tafel  slopes. 

1 .  Four  Point  Central  Difference  Method 

When  the  data  paints  are  equally  spaced  the  Central 
difference  method  may  be  used  to  approximate  the  derivatives 
at  a  point.  The  use  of  this  method  can  result  in 
considerable  time  savings  over  the  cubic  spline  method  as 
the  derivative  can  be  calculated  without  having  to  resort  to 


the  salving  of  a  system  of  equations  for  each  point  in 
question.  The  method  is  not  as  accurate  as  the  cubic  spline 
method  and  graphical  representations  of  the  derivatives  will 
often  contain  more  "noise"  than  those  derived  from  the  cubic 
spline  method. 

The  four  point  Central  Difference  formula  for 
computing  derivatives  is: 


_ d^^ _ _  =  _ 12h _ _ _ 

d(logio  ict)  a(  -  f-r)  -  Va) 

Where : 

h  =  0±  —  01-f]. 

f„  =  logxo(ic«)  at  the  I  th  location  from 
1  oq  X  o  (  i  el )  I 

The  Central  Difference  method  has  an  error  term 
associated  with  it.  In  this  case  the  error  term  is  on  the 
order  of  1/h^.  Although  this  may  appear  to  be  a  substantial 
term,  it  should  be  remembered  it  is  a  possible  error,  the 
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method  is  only  being  used  to  determine  the  Tafel  region  and 
when  used  the  Tafel  constant  is  inverted  so  the  error  term 
would  actually  be  h**. 

2 .  Cubic  Spline  Interpolating  Polynomial 

When  trying  to  numerically  determine  derivatives  of 
experimental  data  which  may  or  may  not  be  equally  spaced  the 
use  of  interpolating  polynomials  is  often  suggested.  Qne  of 
the  most  accurate  methods  is  to  pass  a  cubic  spline  thru  the 
data  points  and  then  to  dif ferentiate  the  resulting  equation 
to  arrive  at  the  derivative  at  a  desired  point.  [Ref.  18: 
pp.  242] 

The  advantage  of  using  a  cubic  spline  is  that  in 
deriving  the  system  of  equations  to  arrive  at  the  constants 
a, b, c,  and  d  which  form  the  equation: 

y  =  ai(x  -  xx)’  -»■  bi(x  -  xi)=  +  Ci(x  -  Xx)  +  dx 
the  system  is  formed  in  such  a  manner  that  the  value 
of  the  function,  first  derivative  and  second  derivative  are 
the  same  for  the  pair  of  cubics  which  join  at  each  point. 
This  requires  that: 

y'  =  3ax(x  -  Xx)^  *  2bx(x  -  xx)  +  c, 
and  y'  '  =  ba,  ( x  -  x,)  +  2bx 

In  order  to  simplify  the  the  mathematics  involved 
the  equations  for  the  system  are  written  in  terms  of  the 
second  derivatives  of  the  interpolating  cubic  Si. 
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The  resulting  interpolating  cubic  will  then  take  on 
the  form  of: 

hi  —  iSi  —  i.  +2S I  (  h  1  —  1 +S  X  )  ■*'h  X  S 1  ♦  (yi»i~yx  )~  (yx~"yi  —  i) 

hx  hi— X 

where  hx  =  x  -  Xi 

Solving  the  system  leads  to  the  following  operations 

to  determine  the  values  of  ax , bx , Cx  and  dx  . 

ax  =  (Sx-«-i  ~  Sx)/bhx 
bi  =  Sx/2 

cTi  =  ((yi-».x  ~  yx)/hx)  ~  ((2hxSx  hxSx-^x)/b) 

dx  =  yx 


The  resulting  system  of  n  -  2  equations  in  Si 
involves  n  pairs  of  data  points  in  order  to  generate  the 
required  number  of  equations.  To  arrive  at  the  two 
additional  equations  to  solve  for  Si  and  S„,  constraints  are 
specified  which  pertain  to  the  conditions  at  the  ends  of  the 
curves.  The  three  choices  for  the  end  conditions  are: 


1.  Si  =  S„  =  O,  which  implies  that  the  end  cubics 
approach  linearity  at  their  extremities. 

2.  Si  =  S3,  Sr,  =  Sr,-i.  This  assumes  that  the  end  cubics 
approach  parabolas  at  their  extremities. 

3.  Si  is  a  linear  extrapolation  of  Si  and  Ss-  S„  is  a 
linear  extrapolation  of  Sr,_3  and  S^-i . 
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In  this  case  since  the  interest  is  in  determining 
the  derivatives  d0/d(logxo  ic«)  pointwise,  and  not  in 
developing  an  interpolating  polynomial,  it  Mas  decided  to 
use  the  end  conditions  of  Si.  =  Sn  =  O,  and  to  use  enough 
data  points  as  to  force  the  inconsistences  caused  by  the 
assumed  end  conditions  aMay  from  the  point  in  question. 

Four  data  points  on  each  side  of  the  point  in  question  were 
used  in  determining  the  the  derivative  at  the  point  in 
question.  This  resulted  in  the  solving  of  a  7  x  7  system  of 
equations  for  each  data  point.  By  locating  the  point  in 
question  in  the  center  of  the  data  set,  the  accuracy  of  the 
estimate  involved  Mas  improved,  and  the  mathematics  involved 
in  determining  the  derivative  Mere  simply  solving  for  c±. 

^ •  Linear  Regression 

Once  the  anodic  and  cathodic  Tafel  regions  have  been 
determined,  the  data  points  Mhich  comprise  these  regions  are 
used  in  a  least  squares  curve  fitting  routine  to  determine 
an  equation  of  the  form: 

0  =  ^(logxo  (in))  +  b, 

Mhere  0  represents  the  Tafel  slope.  The  values  of  a 
and  b  are  found  by  solving  the  tMO  folloMing  equations 
simultaneously. 

logio(  in  )  X  )^  +  ££(logio(  in)  I  )  =  E(  logi.o(  in  )  i  )  ( 0x  ) 

PZ  (  1  og  i.o(in)i)  fW  =  Z  ( 0 1  ) 

Mhere  N  is  the  number  of  points  in  the  Tafel  region. 


5. 


Considerations  Affecting  Method  Selection 
When  the  PAR  351  Corrosion  Measurement  system  scans 
the  potential  durinq  a  Poten tiodynamic  test  the  current  is 
measured  at  intervals  of  2mV  or  0.5mV  of  potential, 
depending  on  the  relative  change  in  the  current.  The 
sampling  rate  is  normally  0.5mV  when  the  current  shows  a 
substantial  change  and  2mV  when  it  does  not.  This  usually 
results  in  a  sampling  rate  of  0.5mV/  in  the  Tafel  regions, 
although  severe  concentration  polarization  in  the  vicinity 
of  ie:e>.-i-  may  result  in  a  limited  number  of  points  having 
O.SmV  intervals  on  the  cathodic  branch.  Two  other  aspects 
may  also  effect  the  potential  intervals  at  which  the  current 
appears  to  have  been  sampled. 

When  data  is  sent  out  the  serial  port,  the  351 
normally  uses  a  format  similar  to  the  FORTRAN  Ga.b  type. 

Not  counting  leading  zeros,  four  significant  figures  are 
used.  This  results  in  a  problem  in  the  recorded  potential 
values.  When  the  potential  is  greater  than  or  equal  to  1.0 
volts,  or  less  than  or  equal  to  -1.0  volts,  only  three 
decimal  places  are  carried.  As  a  result,  when  the  values  of 
0c:ar-r-  fall  in  this  range  and  the  scan  rate  is  0.5mV,  the 
recorded  data  seems  to  reflect  that  the  potential  has  not 
changed  for  two  successive  current  measurements. 

Since  the  potential  is  being  scanned  it  is  rather 
obvious  when  this  occurs  that  the  recorded  data  does  not 
accurately  represent  the  enperimental  data.  Since  this 


pattern  is  easily  recognized,  the  array  containing  the  data 
can  be  manipulated  to  more  accurately  reflect  the 
experimental  data. 

Occasionally  the  potentials  recorded  display  no 
regular  pattern  and  the  data,  although  questionable  in 
nature,  must  be  used  as  recorded. 

When  calculating  the  derivatives,  the  anodic  and 
cathodic  branches  are  performed  separately.  The 
calculations  begin  at  the  third  data  point  aMay  from  Eeot->- 
Mhen  the  Central  Difference  method  is  used,  and  at  the  fifth 
Mhen  using  the  cubic  spline  method.  The  calculations 
continue  as  long  as  the  data  points  continue  to  suggest  a 
slope  of  of  the  expected  sign  (+/-).  This  automatically 
excludes  the  passivating  region  of  the  anodic  section,  as  xt 
should  since  it  is  not  part  of  the  Tafel  region,  and 
excludes  the  possibility  of  the  program  terminating 
execution  if  it  were  to  calculate  an  infinite  slope.  In 
general  this  strategy  worked  as  expected  and  provided  enough 
data  points  to  properly  examine  the  Tafel  regions.  The  only 
experiment  in  which  it  failed  was  for  the  50N05T0N  sample 
tested  in  3.57.  NaCl  solution.  As  can  be  seen  by  examining 
the  Potentiodynamic  curve  of  this  sample  (Figure  3),  the 
cathodic  branch  exhibited  behavior  throughout  which  caused 
either  method  to  terminate  execution  without  generating  any 
derivatives.  The  reason  for  this  behavior  is  unknown. 


The  program  Mas  developed  so  as  to  allow  for  its  use 
without  the  supposition  that  the  user  was  familar  with  the 
FORTRAN  langauge  or  with  personal  computers  in  general. 
Numerous  efforts  were  made  to  minimze  the  possibilty  of  the 
program  prematurely  terminating  execution  in  the  event  of 
improper  input  or  through  floating-point  operations 
resulting  in  invalid  expressions.  By  takxng  advantage  of 
the  color  capability  of  the  monitor  screen,  colors  are 
changed  with  respect  to  the  type  of  input  expected,  and  to 
display  messages  resulting  from  either  input  errors,  data 
file  format  errors  or  a  failure  to  properly  follow  the  steps 
necessary  for  proper  program  execution. 

The  program  was  written  using  Microsoft  FORTRAN 
Optimizing  Compiler  version  4.0.  Graphics  were  achieved 
using  Microcompatibles,  Inc.'s  GRAFMATIC  library.  Plotter 
support  was  achieved  using  Microcompatibles,  Inc.'s 
PLOTMATIC  library. 

A  summary  of  the  main  program  and  its  subroutines  is 
given  below.  The  main  program  and  its  called  subroutines 
are  listed  in  Appendix  C.  The  subroutines  are  presented  in 
the  order  in  which  they  appear  in  the  main  program. 


A .  START . FOR 


This  is  the  main  program.  The  name  of  the  first  data 
file  to  be  used  is  entered.  Following  calls  to  OATAIN, 

CHECK  and  DATAM  a  menu  consisting  of  6  options  is  presented: 

1.  Plot  the  potentiodynamic  curve 

2.  Generate  Tafel  Slopes  using  cubic  spline 

3.  Generate  Tafel  slopes  using  central  difference 

4.  Use  another  data  file 

5.  Overlay  two  potiodynamic  curves 

6.  Exit 

Option  1  generates  calls  to  GRAPHl  and  PLOT 

Option  2  calls  SLOPE 1 .  Following  this  a  first  order 
linear  regression  of  the  anodic  and  cathodic  branches  using 
locations  within  the  original  array  determined  during  a  call 
to  GRAPHl  is  performed  and  the  results  along  with  original 
potentiodynamic  curve  are  passed  to  GRAPHl. 

Option  3  calls  SLOPE.  Upon  return  from  SLOPE  this 
option  proceeds  as  Option  2 

Option  4  returns  the  user  to  the  3rd  executable  line  of 
START. 

Option  5  calls  DATAIN,  CHECK  and  DATAM  while  retaining 
the  data  from  the  first  file. 

Option  6  obviously  exits  the  program. 

B.  DATAIN.FOR 

This  opens  the  input  file  and,  using  a  formatted  read 
statement,  generates  the  array  containing  potential  and 
current  density.  Additionally  this  subroutine  counts  the 
total  number  of  lines  in  the  file.  This  is  done  to  ensure 


that  formatting  changes  within  the  data  file  which  were  not 


corrected  during  the  data  file  transfer  will  not  have  not 
caused  the  array  generation  to  prematurely  terminate.  If  a 
formatting  error  has  been  detected,  an  error  message  is 
displayed  in  red  and  the  line  number  within  the  file  is 
identified.  This  allows  the  user  to  exit  the  program  and 
correct  the  error  before  continuing. 

C.  DAT AM. FOR 

This  converts  the  absolute  value  of  the  current  density 
to  its  logxo  value  and  determines  the  position  of  with 

the  array. 

D.  CHECK. FOR 

This  subroutine  determines  the  number  of  data  points  in 
the  anodic  and  cathodic  branches  to  be  used  in  the 
determining  the  derivatives  when  using  the  central 
difference  method.  This  is  done  by  checking  log  current 
density  for  values  which  might  result  in  the  calculation  of 
an  infinite  slope  (normally  the  onset  of  pitting  or 
passivation)  and  for  equal  increments  in  the  recorded 
potentials.  Occasionally  the  recorded  potentials  will  not 
reflect  the  polarization  which  has  taken  place.  This  occurs 
when  the  absolute  value  of  the  recorded  potential  is  equal 
to  or  greater  than  1.0.  In  these  ranges  the  0.5mV  voltage 
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increment  falls  outside  of  the  total  Midth  of  the  field  and 
is  not  recorded.  As  this  is  a  recognizable  pattern,  at  the 
users  option  the  problem  can  be  corrected. 


E .  SLOPE . FOR 

This  is  the  first  subroutine  called  when  using  the 
Central  Difference  method.  It  assigns  graph  titles,  axes 
labels  and  calls  DAT ADEL. FOR. 


F .  DATADEL . FOR 

Using  the  points  determined  in  CHECK,  this  calculates 
the  numerical  derivatives  using  the  Central  Difference 
method . 


G.  SLOPE 1. FOR 

This  assigns  graph  titles  and  labels.  It  also 
determines  the  number  of  points,  HI,  to  be  used  on  both 
branches  in  the  cubic  spline  method  by  checking  only  for 
values  of  log^o  current  density  Mhich  could  result  in  the 
calculation  of  slopes  of  infinite  value.  It  calls  CSPLIN. 
separately  for  the  anodic  and  cathodic  branches.  MI  anodic 
is  not  necessarily  egual  to  MI  cathodic. 


H.  CSPLIN. FOR 

For  MI  number  of  points,  this  generates  the  7x7  matrix 
and  the  7x1  vector  of  values  described  in  the  cubic  spline 
section  MI -9  times.  As  each  matrix  is  formed  it  is  solved 
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by  a  call  to  LINSYl.FOR.  Upon  return  from  LINSYl  the  7x1 
vector  contains  the  solution  which  is  used  to  calculate  the 
derivative. 


I.  LINSYl.FOR 

Using  partial  pivoting  this  subroutine  solves  the 
augmented  7x8  array  by  Gaussian  Elimination.  The  original 
7x1  vector  contains  the  solution  CRef.  19]. 


J.  GRAPHl.FOR 

This  is  the  graphics  output  of  the  selected  main  menu 
option.  Through  the  use  of  the  numeric  keypad  as  softkeys 
the  user  can  scale  the  axes  as  desired.  The  option  to 
return  to  the  original  graph  or  the  previous  display  always 
exists.  When  plotting  the  results  of  either  numeric 
di f ferentation  method,  only  those  points  which  consecutively 
fall  within  the  desired  region  are  plotted.  Once  the  axes 
have  been  scaled  such  that  graph  contains  only  the  linear 
region  of  the  anodic  or  cathodic  branches  the  position  of 
the  start  and  stop  points  within  the  original  array 
containing  the  derivatives  are  saved.  Since,  when 
calculating  the  derivatives,  the  data  point  containing  Ecorr 
is  always  used  as  a  reference,  ie.  (depending  on  the  method) 
the  first  derivative  for  each  initial  branch  is  either  3  or 
5  points  away  from  Ecoi'.-,  the  start  and  stop  points  also 
represent  the  position  of  the  Tafel  regions  within  the 


original  potential  and  loglO  current  density  array.  When 


the  graph  is  not  of  the  numerical  derivatives,  any  data 
point  which  falls  within  the  desired  axes  is  platted.  If 
desired  GRAPHl  calls  PLOT. FOR  for  a  hard  copy  output  ot  what 
is  currently  displayed  on  the  screen. 

K.  PLOT. FOR 

This  plots  the  graph  displayed  on  the  screen.  Although 
many  of  the  arguments  passed  from  GRAPHl  are  used  directly 
as  passed,  this  is  not  a  'screen  dump' .  The  resolution  of 
the  plotter  is  much  higher  than  that  of  the  screen.  As  a 
result  what  appears  to  be  a  one  pixel  Jump  on  the  screen 
will  not  plot  as  such.  Two  sets  of  tabulated  data  may  be 
plotted  on  each  set  of  axes.  At  the  users  discretion 
different  color  pens  may  be  used  to  plot  the  axes,  first 
output,  second  output,  and  labels. 

L.  CONFIGURATION 

In  its  final  form,  the  program  has  been  compilied  and 
linked  into  an  executable  file  which  runs  on  a  Zenith  248 
computer  which  is  equipped  with  a  80286  microprocessor,  an 
80287  Numeric  Processor  Extension,  a  EGA  graphics  card  with 
64K  video  RAN,  and  a  color  monitor  which  supports  the  EGA 
color  mode.  The  plotter  used  is  an  Houston  Instruments  DHP- 
40  series  digital  plotter. 

The  plotter  library  is  plotter  series  specific.  As  a 
result  of  this,  the  use  of  a  plotter  other  than  a  Houston 
Instrument's  DI1P-29  or  greater  would  require  obtaining 


another  library  from  liicrocompatibles  Inc.  As  many  of  the 
subroutines  contained  Mithin  PLOTIiATICS  are  similar  to  those 
called  Mithin  GRAFMATICS,  it  can  probably  be  safely  assumed 
that  little  if  any  editing  of  the  subroutine  PLOT. FOR  would 
be  required. 

To  use  the  program  with  a  video  configuration  other  than 
that  previously  described  would  require  editing  of  the 
subroutines  PLOT. FOR,  GRAPHl.FOR,  DATAIN.FOR,  ST ART. FOR  and 
CHECK. FOR.  DATAIN,  START  and  CHECK  would  require  editing 
only  if  the  video  configuration  was  monochrome.  PLOT  and 
GRAPHl  would  require  editing  for  any  graphics  configuration 


other  than  EGA. 


V .  RESULTS  AND  DISCUSSION 

A.  SYNTHETIC  SEAWATER 

Table  1  below  compares  the  Tafel  constants  determined  by 
using  the  new  program  with  those  obtained  by  Escue  [Ref.  9]. 

TABLE  1 

TAFEL  CONSTANTS  OBTAINED  FROM 
SYNTHETIC  SEAWATER  EXPOSURE 


MATERIAL 

0. 

0-* 

0.= 

0=- 

Ti-50Ni 

0.1400 

0.174 

0.1155 

0.132 

1020  Steel 

0.0758 

0.068 

0.1481 

0.114 

304  Steel 

0.0241 

0.178 

0.0207 

7075  A1 

0 . 3790 

0.087 

0.1584 

0.062 

Fe-Cr-Mo 

0.1775 

0.121 

0.1226 

0.141 

Fe-Cr-Al 

0.1840 

0.180 

0.1369 

0.151 

Cu-Mn-Al-Fe-Ni 

0.0399 

0.048 

0.1579 

0.120 

Cu-Mn-Al 

0.0633 

0.035 

0.0652 

0.045 

Cu-Zn-Al 

0.0804 

0.062 

0.0785 

0.075 

630  Bronze 

0.0403 

0.031 

0.0724 

0.041 

t  From  Escue  [Ref.  9] 

Table  2  below  compares  the  current  densities  and 
related  corrosion  rates  of  the  samples  tested  in  synthetic 
seawater.  The  first  line  for  each  material  represents  the 
current  density  as  determined  by  the  intersection  of  the 
Tafel  lines.  In  many  of  the  cases  the  anodic  and  cathodic 
Tafel  lines  did  not  intersect  at  a  single  point  and  are 


.V, 


'••.V 

V  V  V  v'  N’-V' 


joined  by  a  horizontal  line.  In  these  cases  ieoi-r  was 
determined  by  taking  the  average  of  the  anodic  and  cathodic 
icc3i-r-  at  the  point  where  the  lines  intersected  a  horizontal 
line  form  by  deor-i-. 

The  second  line  represents  the  results  obtained  by 
combining  the  Tafel  slopes  with  the  polarization  resistance 
data. 

TABLE  2 


CURRENT  DENSITY  AND  CORROSION  RATES 


MATERIAL 

ico.'.' 

ic  or-r- 

CR 

CR» 

uA/cm* 

uA/cm^ 

mpy 

mpy 

Ti  -  SO*/.Ni 

1.245 

1.734 

0.573 

0.799 

4.182 

4.960 

1.822 

2.160 

1020  Steel 

4.212 

3.837 

1.981 

1.805 

20.325 

17.201 

9.580 

8.110 

7075  A1 

1.413 

0.421 

0.591 

0.176 

3.067 

0.994 

1.281 

0.415 

304  Steel 

0.106 

0.728 

0.048 

0.329 

0.157 

1.190 

0.071 

0 . 537 

Fe-Cr — Mo 

0.750 

0.804 

0.354 

0.380 

2.294 

2.060 

1.085 

0.974 

Fe-Cr-Al 

0.822 

0.830 

0.362 

0.365 

1.183 

1.240 

0.520 

0.545 

Cu— Mn- A 1 — Fe— 

Ni  8.241 

8.854 

4.165 

4.070 

9.662 

10.400 

4.882 

5.260 

Cu-Mn-Al 

2.089 

1.119 

1.111 

0.595 

1.459 

0.894 

0.784 

0,481c 

Cu-Zn-Al 

2.239 

1.897 

1.136 

0.963 

3.780 

3.230 

1.919 

1.640 

630  Bronze 

1.762 

1.442 

0.901 

0.737 

•  From  Escue 

6.878 
CRef.  91 

4.690 

3.550 

2.420 

Ti-50y.Ni  (TINILb) 


These  results  compared  favorably  with  those  obtained 


by  Escue.  Figure  4  shows  the  derivatives  as  calculated  by 


the  cubic  spline  method.  As  can  be  seen  in  this  figure  no 


true  linear  regions  are  distinguishable  in  either  of  the  two 


branches.  When  this  occurs  the  "knee"  method  can  be 


employed  by  examining  the  branches  for  a  range  where  the 


derivatives  oscillate  about  a  linear  midpoint.  In  this  case 


bot  '  he  anodic  and  cathodic  branches  showed  this 


oscillation  in  the  log  current  density  range  of  —5.28  to  - 


5.51.  Figure  5  shows  the  resultant  Tafel  lines  plotted  on 


the  original  potentiodynamic  curve.  Although  the  curves 


intersect  dec.-,-  at  what  could  be  considered  to  be  a  single 


point  the  icor-t-  determined  in  this  manner  differs  greatly 


from  that  calculated  by  using  the  Tafel  constants  and  the 


polarization  resistance  data. 


2 .  1 PJZOL  Ca r bon  .St ee L.. ( S1.P2PLG ) 


Figures  6  and  7  show  the  results  as  obtained  for  the 


1020  carbon  steel  sample.  Of  all  of  the  samples  tested  this 


one  exhibited  a  cathodic  branch  with  the  least  tendency  to 


behave  in  a  linear  fashion.  If  tnere  was  any  linear 


behavior  in  the  cathodic  branch  it  occured  at  a  logi.o  ico.-.- 


of  about  -5.2.  The  anodic  branch  displayed  a  much  longer 


lasting  region  around  logi.o  ic 


=  -  4.3.  This  sample  also 


showed  the  greatest  disparity  in  the  value  of  i^cDf-r 


depending  on  the  manner  in  which  it  was  calculated. 


r,  v"-  ►"v  •' .  ' 


3 .  7P?5  Aluminun,  ( AL7p7.5L6 ) 

In  this  case  the  central  difference  method  showed  a 
linear  region  in  both  the  anodic  and  cathodic  branches.  As 
can  be  seen  in  Figure  0  the  linear  regions  exist  at  the 
pronounced  peaks  in  the  derivatives.  In  this  case  the  Tafel 
lines  did  intersect  at  a  single  point  as  seen  in  Figure  9. 

4 .  304  Stainless  Steel  (SST304LG) 

The  calculated  corrosion  rates  and  the  respective 
Tafel  constants  differed  greatly  from  those  obtained  by 
Escue.  Although  the  calculated  corrosion  rates  correlated 
better  to  the  actual  weight  loss  methods,  the  need  to 
terminate  the  calculating  of  derivatives  before  the  onset  of 
passivity  may  have  caused  a  linear  region  in  the  anodic 
branch  to  be  omitted.  These  points  can  be  seen  by  the  sharp 
peak  on  the  anodic  branch  on  Figure  10,  and  the 
corresponding  point  on  the  poten tiodynamic  curve  Figure  11. 
The  Tafel  lines  intersected  at  a  single  value  of  log^o  icroi-r- 
but  this  time  there  was  a  good  correlation  between  the 
calculated  and  graphical  icor-r-* 

5  •  FeHCr -Mq_ ( VACRO ILGl 

Again  the  presence  of  any  truly  linear  regions  is 
questionable,  although  both  branches.  Figure  12,  show  a 
range  where  the  deviations  from  linear  behavior  are  minimal. 
Figure  13  once  again  shows  a  theoretical  intersection  of  the 
Tafel  lines. 


6. 


Ft?-Cr-Al  (yACRQ2LG) 


This  alloy  showed  behavior  similar  to  the  Fe-Cr — Mo 
alloy  described  above.  What  may  be  considered  a  truly 
linear  region  in  the  cathodic  branch  is  shown  in  Figure  14 
just  to  the  left  of  the  sharp  peak.  The  overlay  of  the 
Tafel  lines  on  the  potentiodynamic  curve  is  shown  in  Figure 
15.  The  Tafel  constants  derived  show  a  closer  relationship 
ho  those  calculated  by  Escue,  than  did  a  comparision  of  the 
Fe— Cr-Mo  sample. 

7  -  «n -A 1 rfer  N i  ( saNOSTLG ) 

SQNQSTQN,  Figures  16  and  17,  displayed  a  linear 
anodic  branch  and,  as  did  most  of  the  alloys,  a  cathodic 
branch  effected  by  concentration  polarization. 

8.  Cu-Mn-Al  (INCRLG) 

The  result  here  again  showed  a  linear  anodic  branch 
but  in  this  case  the  concentration  polarization  did  not 
overwhelm  the  cathodic  branch  in  such  a  manner  as  to  distort 
the  region  where  activation  polarization  is  predominant. 
Figure  18  shows  the  graph  of  all  the  derivatives  and  Figures 
19  and  20  show  the  graphs  of  the  derivatives  in  what  was 
decided  to  be  the  Tafel  region.  Figure  21  shows  the 
original  potentiodynamic  and  the  Tafel  lines. 

9 .  Cu- Z  n - A  1  ( DL ALCLG ) 

Experinced  similar  behavior  as  the  Cu— Mn— A1  sample. 
Once  again  a  linear  region  was  very  apparent  in  the  anodic 
branch.  Figures  22  and  23. 


10.  630  Series  Bronze  (BRNZLG) 

This  Mas  the  only  sample  of  those  tested  in  the 
synthetic  seaMater  where  the  Tafel  reqions  for  the  anodic 
and  cathodic  branches  were  found  in  the  same  region  of  log 
ieof-r-.  As  seen  in  Figure  24  concentration  polarization 
again  quickly  dominates  the  cathodic  branch.  Figure  25 
displays  the  overlay. 

B.  3.57.  NaCl  SOLUTION 

In  a  similar  fashion  Table  3  lists  the  Tafel  constants 
and  Table  4  the  resultant  current  densities  obtained  from 
the  samples  tested  in  the  3.5iC  NaCl  solution.  In  general 
the  current  results  show  a  stronger  correlation  with  the 
corrosion  rates  obtained  by  Akthar  [Ref.  lO]  than  did  the 
comparision  of  the  current  results  with  those  obtained  by 
Escue  [Ref.  9]  for  the  synthetic  seawater. 

TABLE  3 

TAFEL  CONSTANTS  OBTAINED  FROM 
3.5  7.  NaCl  SOLUTION  EXPOSURE 


MATERIAL 

B. 

Be 

Be* 

304  Steel 

0.1259 

0.4117 

0.0951 

0.117 

7075  A1 

0.0116 

0.010 

0.0303 

2.36 

Fe— Cr — Mo 

0.2821 

0.234 

0.1053 

0.110 

Fe-Cr-Al 

0.1203 

0.214 

0.1292 

0.1205 

Cu-Mn-Al 

O . 0002 

0.0137 

0.0337 

1.379 

630  Bronze 
*  From  Akthar 

0.0524 
[Ref.  103 

0.570 

0.1472 

0.321 

TABLE  4 


CURRENT  DENSITY  AND  CORROSION  RATES 


MATERIAL 

A  e:c^r“r- 

« 

L 

L 

0 

U 

■H 

CR 

CR* 

uA/cm= 

uA/cm=* 

mpy 

mpy 

7075  A1 

6.960 

7.406 

0.048 

2.965 

10.460 

11.700 

1.842 

4.685 

304  Steel 

0.052 

0.100 

0.026 

0.050 

0.012 

0.020 

0.006 

0.010 

Fe-Cr — Mo 

0.043 

0.040 

0.019 

0.018 

0.609 

0.613 

0.274 

0.276 

Fe-Cr-Al 

0.493 

0.660 

0.234 

0.313 

0  ‘'■.6 

1.180 

0.428 

0.551 

Cu— Mn— A1 

2.864 

6.062 

1.454 

3.078 

5.251 

10.800 

2.666 

5.483 

630  Bronze 

0.093 

0.129 

0.048 

0.066 

3.570 

4.480 

1.842 

2.312 

•Fro*  Akthar  QRef.  10] 

1 •  7075  Aluminum  (AL7075) 


In  this  case  the  effect  of  concentration 
polarizaticsn  made  any  estimate  of  the  cathodic  Tafel 
constant  extremely  difficult.  As  seen  in  Figures  26  and  27 
concentration  polarization  completely  overwhelmed  the 
cathodic  branch.  As  a  result  only  the  "knee”  method  cuold 
be  used,  but  its  accuarcy  is  questionable.  The  anodic 
branch  displayed  typical  linear  behavior. 

2 -  504  Stainless  Steel  (SST504 ) 

The  cathodic  branch  of  this  sample  did  display  a 
linear  region  located  between  the  sharp  peaks  in  Figure  28. 
The  anodic  branch  displayed  numerous  attempts  to  passivify 
as  seen  in  Figure  29. 


i  3 .  Fe-Cr-Mo  (VCMQPD) 

I 


Similar  behavior  to  the  304  Stainless  steel  above. 
The  bandwidth  exhibited  by  the  cathodic  branch  near  logxo 
ic:ci.-t-  of  about  —  6.5  Figure  30,  typifies  the  knee  in  the 
curve.  Again  numerous  attempts  to  passivity  limited  the 
range  of  data  which  could  be  use  in  the  anodic  branch  Figure 


31. 


^ •  Fe-Cr-Al  (VACAL2) 

Both  the  anodic  and  cathodic  regions.  Figure  32, 
displayed  a  region  where  the  behavior  could  be  considered  to 
be  linear.  Although  the  Fe— CR— Mo  and  Fe— Cr — A1  samples 
tested  in  synthetic  seawater  displayed  similar  behavior,  the 
Fe-Cr — A1  sample  tested  in  the  3.5%  NaCl  solution  did  not 
exhibit  the  tendency  to  passivity.  Figure  33,  as  did  the  Fe- 
Cr-Mo  sample. 

5 •  Cu-Mn-Al  (INCRMTE) 

This  sample  best  typifies  the  problem  associated 
with  concentration  polarization.  On  Figure  34  a  linear 
region  can  be  seen  in  the  cathodic  branch  at  a  logxo  ieoi-r 
of  about  -5.3.  This  region  is  followed  by  a  rapid  decrease 
in  the  slope  of  the  curve  indicative  of  concentration 
polarization.  When  the  original  potentiodynamic  curve  is 
viewed  thru  out  its  entire  polarization  range.  Figure  35, 
very  little  can  be  determined  due  to  the  great  changes  in 


current  density  on  the  anodic  branch.  When  the  range  in  the 


cathodic  branch  displaying  the  linear  bevaior  is  scaled 


I 


properly  a  textbook  example  of  concentration  polarization  is 
seen.  Figure  36.  The  overlay  is  shown  in  Figure  37. 

6 .  630  Series  Bronze  (BRZ630PD) 

As  seen  on  Figures  38  and  39  this  sample  experienced 
similar  behavior  to  the  one  above.  All  though  the  linear 
region  directly  before  concentration  polarization  is  not  as 
well  defined  the  cathodic  behavior  is  dominated  first  by 
activation  polarization  and  then  by  concentration 
polarization . 


C.  COnPARISION  OF  CURRENT  DENSITIES 

In  many  of  the  samples  tested  the  corrosion  rates  as 
determined  by  the  intersection  of  the  Tafel  lines  were 
significantly  lower  than  those  obtained  by  the  use  of  the 
Tafel  Slopes  and  linear  polarization  data.  As  previously 
stated  the  Tafel  regions  should  begin  to  dominant  the 
polarization  curve  at  an  polarization  increment  of  about 
50mV  from  Ecot-f  This  implies  that  a  linear  region  should 
not  be  present  at  increments  less  than  this.  With  these 
tests  a  linear  region  was  often  found  within  this  zone. 
Because  of  the  general  shape  of  a  potent iodynamic 
polarization  curve  the  resultant  lines  drawn  from  this 
region  will  always  have  a  lower  slope  and  thus  intersect  at 
a  lower  value  of  i.  The  disparity  underscores  the  need  to 
use  the  average  of  as  determined  by  the  intersection 


of  the  Tafel  lines  and  that  calculated  by  the  Tafel  slopes 
and  linear  polarization  data  when  estimating  corrosion 
rates . 

Another  interesting  feature  was  noticed  when  conducting 
some  data  manipulation  of  the  potentiodynamic  polarizations 
analyzed.  When  the  absolute  value  of  (0  -  0c:e>'r-)  was 
plotted  against  log^o  ia  in  many  of  the  cases  the  anodic  and 
cathodic  branches  began  to  deviate  significantly  at  a  log  ia 
corresponding  to  the  intersection  of  the  two  Tafel  lines. 
This  may  be  coinicidental  but  also  may  warrant  further 
discussion.  When  a  linear  polarization  resistance  test  is 
conducted  the  potential  is  scanned  with  a  range  of 
25  mV.  The  polarization  resistance  is  determined  by 
calculating  the  slope  of  the  linear  region  about  i  =  O 
A/cm2.  If  a  linear  polarization  curve  is  plotted  in  the 
same  manner  as  above,  while  neglecting  the  singularity  at  i 
==  O  A/cm^,  the  deviation  discussed  would  determine  the  end 
of  the  linear  region.  Returning  to  the  equation  used  to 
determine  iee>.-r  using  the  Tafel  constants  and  the 
polarization  resistance  it  can  be  shown  that  in  order  to  the 
polarization  resistance  to  be  a  constant  that  PsPc  = 

where  C  is  a  constant  of  proportionality.  The 
current  being  measured  during  a  potentiodynamic  polarization 
IS  equal  to  the  difference  between  the  anodic  and  cathodic 
currents.  Only  when  one  of  the  currents  begins  to  overwhelm 
the  other  would  there  be  a  substantial  change  in  the 


measured  current  Mith  respect  to  a  given  overvoltage.  This 
would  correlate  to  the  deviation  between  the  anodic  and 
cathodic  branches  on  a  plot  of  the  absolute  value  of  (0  - 
0cor>-)  vs  log  (id).  It  also  should  correspond  tc  the  end 
of  the  linear  region  on  the  linear  polarization  curve  since 
at  this  point  the  constant  of  proportionality  wcjuld  no 
longer  be  a  constant  as  one  of  the  two  branches  is  beginning 
to  overwhelm  the  other.  It  is  postulated  that  a  technique 
of  this  method  could  be  used  to  determine  the  minumium 
overvoltage  at  which  the  search  for  a  Tafel  constant  could 
begin.  It  is  recognized  that  this  does  not  account  for  the 
possibility  of  completely  symetric  anodic  and  cathodic 
branches  but  if  this  was  the  case  there  would  be  no  problem 
in  resorting  to  the  +/-  50mV  rule  proposed  by  Ailor. 


VI,  COf^CLUSigN  AND  RECOMMENDATIONS 

The  FORTRAN  based  graphics  program  proved  to  be 
extremely  valuable  in  determining  the  existence  of  linear 
Tafel  regions  particularly  in  those  cases  where 
concentration  polarization  dominated  the  cathodic  branch. 

With  the  ability  to  transfer  the  experimental  data  to  a 
programable  computer  the  opportunity  exists  to  further 
analyze  the  data.  This  could  allow  for  the  determination 
which  electrodes  are  effecting  the  main  redox  reaction  when 
the  environment  consists  of  a  solution  of  mixed  electrodes. 
In  order  to  accomplish  this  it  is  recommended  that  a  single 
high  damping  alloy  be  subject  to  individual  poten tiodynamic 
polarizations  containing  a  single  component,  at  its  normal 
concentration,  of  synthetic  seawater.  Through  the  use  of 
graphical  techniques  it  should  be  possible  to  isolate  the 
influential  redox  reactions. 

In  order  to  accomplish  this  it  will  be  necessary  to 
conduct  all  the  experiments  in  solutions  which  contain  the 
same  amount,  or  at  least  a  monitored  amount,  of  dissolved 
oxygen.  This  is  one  factor  which  was  not  considered  in  the 
previous  two  experiments  and  may  have  made  the  analysis  more 
difficult  than  necessary. 
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APPENDIX  B 

PAR  351  TO  2-248  DATA  TRANSFER  PROGRAM 

10  J  =  O 

20  SCREEN  O 

30  COLOR  4,1,1 

40  KEY  OFF;CLS:CLOSE 

50  DEFIMT  A-Z 

60  LOCATE  25,1 

70  PRINT  STRING*(60,"  ") 

80  FALSE  =  0:7RLE=N0r  FALSE 
90  X0FF*=CKf»(19):X0r«==a-R*(17) 

lOO  LOCATE  25,20:PRINr  “FW  351  TD  Z-248  DATA  TRANffER  PROGRAM" 
no  LOCATE  3,20:PRINr  "YOU  MLET  FRDFEPLY  CDNFIOJR  TVt  PAR  351." 

120  LOCATE  5,20:FRINT  "AFTER  SETTING  TFE  Tir*"  TvC  DATE  ON  TFE  351" 

130  LOCATE  7,20:PRINr  "FREBS  SYSTEM  MANAGER  " 

140  LOCATE  8,20:KRINT  "PFESS  CDKFIGLJR  SYSTEM  " 

150  LOCATE  9,20:miNr  "SET  FROTOCOL'  TO  PRINT" 

160  LOCATE  10,20;PRINT  "^T  RARITY'  TD  NOKE" 

170  LOCATE  11,20:FRINT  "SET  BALD  RATE'  TD  2400" 

180  LOCATE  12,20:FRINr  "SET  'STOP  BIT'  TD  1" 

190  LOCATE  13,20:PRINT  “SET  'WORD  LEMTIR'  TO  8" 

•200  COLOR  14,1,1 

210  LOCATE  15,20:PRINT  "PFESS  ANY  KEY  TO  CONTINUE" 

220  GD*  =  IM<EY*:IF  GO»=  ""  GOTO  220 
230  CLS 

240  CCLER  4,1,1 

250  LOCATE  25,20:FRINT  "PAR  351  TO  Z-248  DATA  TRANSFER  PROGRAM" 

260  LOCATE  3,20:PR1NT  "PFESS  'MAIN  fEMJ'  ON  TFE  351" 

270  LOCATE  4,20:PRINT  "PFESS  'RECALL/DISFLAV  EXPERIMENT  " 

280  LOCATE  5,20:PRINT"PFESS  'DI34_AY  EXPERIFCNT'CN  TFE  CDFY  FF01  DISK 
SECTION" 

290  LOCATE  6,20;PRINT  "SELECT  TFE  EXFERIfENT  TO  BE  TFWN3ERED" 

300  LOOTTE  7,20:FRINT  "PRESS  ' PLOT/FCFMAT '  " 

310  LDCAIE  8,20:PRINT  "EIELFE  TFE  X  AXIS  IS  IN  TFE  LIKEPR  SCALE" 

320  COLOR  2,1,1 

330  LOCATE  9,20:PRINT  "IF  NJT  IN  A  LIIEAR  DISPLAY  TFEN  PFESS" 

340  LOCATE  10,20:PRINT  "  FTFMAT  DISPLAY,  LIKEAR' , 'FEVIEW  DISPLAY ',’ RJD 
FORMAT  " 

350  COLOR  14,1,1 

360  LOCATE  14,20:PRINT  "PFESS  ANY  KEY  TD  CONIINE" 

370  GD*=IN<EY*:  IF  QD*=  ""  GOTO  370 
38C'  CLS 

390  COLOR  4,1,1 

400  UXIATE  25,20:FRINT  "PAR  351  TO  Z-248  DATA  7TTANBFER  FRCOVM" 

410  SEED*  =  "2400" 

420  COrt-IL*  ="C0M2:"+SPEED*+",N,8,1,CS,CS" 

430  OFEN  COTEIL*  AS  #1 


I 


4«)  PRINT  1»1,  XCFF*; 

450  CPtN  "SCRN;"  FOR  OJIFUT  AS  4»2 

LXCATC  10,25:PRINT  ■•(R)eceive  a  file  OR  (E)xit  ” 

470  LOCATE  11,37;  I^RJT  TXRXS 

480  IF  (TXRX*<>"R")  AND  ( TXRX«>*‘E" )  TT€N  460 
490  IF  TXRX«=‘*E"  TVEN  1350 
500  CLS 

510  COLOR  4,1,1 

520  LCO«E  25,20;PRINT  "PAR  351  7D  Z-248  DATA  TRANSFER  PROGRAM" 

530  COLOR  15,1,1 

540  LOCATE  13,25:PRINT  •‘PFESB  'PRINT  DATA'  ON  TTE  351  TVEN" 

550  LOCATE  15,30;PRINT  "PRESS  ANT  KEY  TO  CDNTIN£" 

560  GO*  =  INKEY*; IF  G»="“  GOTO  560 
570  CLS 

580  COLOR  4,1,1 

590  LOCATE  25,20;miNT  "PAR  351  TO  Z-248  DATA  TRANSFER  PROGRAM" 

600  LOCATE  10,20;FRINT  "THIS  PHDQTAM  USES  A  RAM  DISK  TO  SPEED 
EXECLITION. " 

610  COLOR  2,1,1 

620  LOCATE  13,20:PR1NT  "DO  NOT  ENTER  TVE  DISK  DRIVE  OR  THE  FILE  TYPE" 
630  LOCATE  14,20;miNT  "BLJT  OM.Y  TVE  NA^E  OF  TVE  OLfTPUT  FILE" 

640  CDLCR  14,1,1 

650  LOCATE  16,20;PRINT  "ENTER  TVE  OLDRIT  FILE  NA»E" 

660  LOCATE  18,25;  IM=LIT  FIL* 

670  D3<*="D;";D6KFIL*=OSK*^IL* 

680  OPEN  "D;TVP"  FOR  OLfTPUT  AS  #3 
690  N  =  0 
700  CLS 

710  PRINT  #1,  XON*; 

720  IF  L0C(l)=O  TVEN  BDSLJB  820 

730  IF  L0C(1)>12B  TVEN  RALBE=TFl£;PRINT  #1,X0FF* 

740  A»=INVr*(l,#l) 

750  IF  (A»=C>«*(10))  OR  (A*=CVF*(13) )  TVEN  B*=CVR*(  13)  :J=J+1 
760  IF  (A*<>OE»(13))  AN)  (A*<XVR*(10) )  TVEN  PRINT  #3,Al*; 

770  IF  (J  >  3)  AND  (B*=CVR»(13) )  AN)  (A*=OR*(13))  TVEN  PRINT 
l»3,CVf»(13) 

780  IF  (J  >  3)  AND  (B*=CVfT*(  13) )  AN)  (A(*=C1E*(13) )  TVEN  B*="  ":N=N«-1 
790  IF  LOC(l)X>  TVEN  730 

000  IF  PAUSE  TVEN  RALJ^=FALSE; PRINT  #1,XDN*; 

010  GOTO  720 

020  FUT  1=1  TO  5000 

830  IF  LOCdXX)  TVEN  I  =  9999 

040  VEXT  I 

850  IF  I>9999  TVEN  FETURN 
860  CLOSE  It3;CLS 
870  F*="D;"*A* 

080  OPEN  "D;TVP"  FOR  IVPUT  AS  #4 
890  OFEN  DSKFIL*  FDR  OUTPUT  AS  #3 
900  FOR  I  =  1  TO  N 
910  X*  =  IVPUT*(44,#4) 

920  Y*  =  IVf4JT*(2,#4) 
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V  J'  ■■  •**  J" 


930  0*  =  MIWCXS.IB,*?) 

940  C*  =  mD»(B*,lfl) 

950  D*  =  MID«(B*,2,1) 

960  E*  =  MID*(B*,3,1) 

970  F*  =  MID*(B*,4,1) 

900  G*  =  mD*(B*,5,l) 

990  H*  =  f1I[»(B*,6,l) 

1000  K  = 

1010  J*  =  MID*(B*,8,1) 

1020  K*  =  mD*(B*,9,l) 

1030  FX*  =  "-0.000000" 

1040  IF  (£♦  =  ”E")  THEN  MIDtlX^.ia**?)  =  FX* 

1050  IF  (F*  =  "E")  TFEN  MID*(X«,1B,9)  =  FX* 

1060  IF  (G*=‘*E'‘)  TIEN  MID*(  X*,  18,*?)  =  FX» 

1070  IF  (H*="E'*)  THEN  riID*(X*,  18,9)  =  FX* 

1000  IF  (H*="E")  TIEN  MID*!  X*,  18,9)  =  FX* 

1090  IF  (I*=‘’E")  TTEN  MID*(X*,18,9)=FX* 

IlOO  IF  (J*="G")  TIEN  MID*(X*,’o  ")=FX* 

lllO  IF  (K*=*'E”)  TFEN  MIDt(X*,18,9)=FX* 

1120  F2*  =  ’^-12  ■■ 

1130  A*  =  MID*(X*,42,1) 

1140  A1*=NID*(X*,43,1) 

1150  A2*=MID*(X*,34,1) 

1160  A3*  =  A*^l* 

1170  IF  (A3*  =  ■•12")  AM)  (A2*  =  •*-")  TFEN  MID*(X*,39,5) 
1100  P*  =  mD*(X*,l,l) 

1190  IF  (P*  =  "P")  TFEN  PRINT#3,X* 

1200  lEXT  I 
1210  (XOEE  #3:(XS 
1220  J  =  0 
1230  (XLCF  1,5,4 
1240  BEEP 

1250  LOCATE  12,30:miNr  "****TRRNSFER  (IIELETE«**$" 

1260  LOCATE  13,30:PRINr  "HESS  ANV  KEY  TD  CCNTIME" 
1270  GD*=IM<EY*:1F  QD*=""  GOTO  1270 
1280  PRINT  *1,  XCFF* 

1290  CIOGE  #1 
1300  C1.0EE  W2 
1310  ooa;  #3 
1320  CLQg  #4 
1330  (OLCF  1,1,1 
1340  GOTO  380 
1350  EM) 


APPENDIX  C 
CORROSION  PROGRAM 


THIS  IS  TVE  MAIN  PROS¥»1.  WITHIN  IT  TVE  USER  ENTERS  TTE  FILE 
NATE  UHICH  hE/JJE  WIShES  TO  ANALYZE.  FOLLOWING  THIS  A  fEMJ 
APPEARS  ALLOWING  THE  USER  TO  SELECT  TFE  DESIFED  OPTICN. 

TFE  FOLLOWING  VARIBLES  ARE  DEFINED 

□DATA  :  TFE  2  DirENSICNAL  ARRAY  FDR  TFE  DATA  FILE  CONTAINING 
POTENTIAL  AND  CURRENT  DENSITY 
CDATAI:  AS  fiOJJE  ELIT  ON-Y  USED  WFEN  CKEFEAYING  TWO  CURVES 
FDATA  :  A  TWO  DIFENSIQNAL  ARRAY  OF  CDATA  BLfT  NOW  IN  TERS  OF 
POTENTIAL  AM)  LOBIO  dJRENT  DENSITY 
XDAT  :  A  O^E  DIFENSICivlAL  ARRAY  CONTAINING  TFE  X  VALUES 
CF  TFE  FIRST  CURVE  TO  EE  GRABED  OR  PLOTTED 
YDAT  :  A  OTE  DIFENBIONAL  AfRAY  CONTAINING  TFE  Y  VALUES 
OF  TFE  FIRST  CURVE  TO  BE  GFWFFED  OR  PLOTTED 
XDATl  :  A  OTE  DIFEIEKNAL  ARRAY  CONTAINING  TFE  X  VALUES 
CF  TFE  SECOM)  CURVE  TO  BE  GRAPFED  CR  FLOTTED 
YDATl  :  A  OFE  DIME^BIONAL  ARRAY  (UNTAINING  TFE  Y  VALUES 
OF  IFE  SECOM)  CURVE  TO  BE  GRAPFED  OR  PLOTTED 
XI  ;  THE  SUM  OF  TFE  LOG  CUJFFENT  DENSITIES.  USED  IN  TFE 

LINEAR  FEGRESSICN  FDR  DETEFtllNING  TFE  TAFEL  CONSTANTS 
X2  :  TFE  SLM  CF  TFE  SQUARES  OF  TFE  LOG  CURRENT  DENSITY 

INITIALLY  USED  IN  THE  SAME  MAMER  A  XI  BUT  CNCE  TFE 
TAFEL  CONSTANT  HAS  BEEN  DETERlIfED  USED  A6  TFE  LOG 
CURFENT  DENSITY  VALUE  IN  CREATING  TFE  TAFEL  LIFES 
XY  ;  TFE  SUM  OF  TFE  FRCDUCTTS  OF  LOG  CUfFENT  DENSITY  AM) 
POTENTIAL.  USED  IN  TFE  SAME  MAMER  AS  XI  AND  X2 
Y1  :  TFE  SLM  OF  TFE  POTENTIAL.  UGED  AS  XI,X2,XY 

NS  ;  TFE  NLTBER  OF  POINTS  USED  TO  GEFERATE  TFE  TAFEL  SLOPE 

B2  :  TFE  TAFEL  SLOPE  AS  CALCULATED  IN  TFE  REGRESSION 

B1  :  TFE  INTERCEPT  CF  TFE  TAFEL  LIFE 

M  :  FENJ  RESPCFEE 

I  :  DO  LOOP  dUNTTER 

N  :  TFE  NLfEER  OF  POIN1S  IN  TFE  DATA  FILE  BEING  RET«) 

01  :  TFE  NUMBER  CF  POINTS  IN  TFE  FIRST  DATA  FILE,  UNCHANGING 

N1  :  TFE  NUrCER  OF  POINTS  IN  TFE  FIRST  DATA  FILE.  USED  IN 

TFE  SUBRDUTIFE  SLOPEl  AM)  CPN  EE  FETUFTED  WITH  A 
DIFFERENT  VALUE 

J  :  TFE  lUCATICW  CF  ECXFR  WITHIN  M)ATA 

J1  :  SAM-  AS  J  aJT  PASSED  AS  AN  ARDUEJENT  TO  SLOFEl 

02  :  SATE  AS  01  BUT  USED  CN-Y  WFEN  CKER-AYIFC  DATA  FILES 

J2  ;  SAFE  AS  J  BUT  USED  CM.Y  WEEN  CKER^YING  DATA  FILES 

SI  :  TFE  POBITICN  WITHIN  FDATA  l*EFE  TFE  FEOUIFEMNIS  CF 

TFE  CENTRAL  DIFFEFMU"  MTHDD  FAIL.  CATHODIC  BRAMH 


52  :  SP^E  AS  SI  BUT  FOR  TVE  ANODIC  BRANCH 

53  ;  TVE  TOTAL  NLrCER  OF  PQIMTS  IN  7T-E  CWTHODIC  BRATCH 

WHim  MEET  TTE  FEUJIREMENTS  FDR  USING  TFE  CENim. 
DIFFERENCE  ^E^H□D 

54  :  SAfE  AS  S3  BUT  FCK  TFE  ANODIC  BRATCH 

A1  :  TFE  POSITION  WITHIN  XDATl  AlO  YDATl  WHICH  MARCS  TFE 
:  BEBIIVNIfJB  OF  TFE  TAFEL  FEBION 
A2  :  SAME  AS  Al  BUT  TFE  EM>  OF  TFE  REGION 

A3  :  SERVES  TFE  S»E  PURPOSE  AS  Al  BUT  WITH  RE3EET  TO 

FCATA 

A4  :  SVE  AS  A3  BUT  AS  A2 

Cl  :  SAfE  AS  Al  BUT  WITHIN  XDAT  AlO  YDAT  (CATFCDIC) 

C2  :  ^3^E  AS  Cl  BUT  TFE  EM)  CF  TFE  REGION 

C3  :  SERVES  TFE  SATE  PURPOSE  AS  Cl  BUT  WITH  RESPECT  TO 

rOATA 

C4  ;  SAfE  AS  C3  BUT  AS  C2 

N0  SET  TO  O.  FWSSED  TO  GRAFHl  AND  PLOT  TO  AMDID  TFE  USE 

OF  TFE  EFPTY  XDATl  AND  YDATl  ARRAYS  WFEN  TIE  AfE  lOT 
USED 

NATE  :  TFE  NAIE  CF  TFE  MAIN  DATA  FILE  IN  USE 
NAME2  :  TFE  NAfE  OF  TFE  SBCOM)  FILE  WFEN  0^€RLAYING 
TITIJE2:  TFE  SECOO  HEADING  FOR  TFE  GfTARH  AM)  PLOT.  CHANGES 
DEPEMJENT  CN  TFE  OPTION  SFl  H:iFD 
TITLES;  RESERVED  FOR  FUTLfE  USE  IF  DEBIFED 
D1  ;  CHARACTER  STRING  USED  TO  FORI  TITLE2 
02  :  TFE  CHARCACTER  REPRESENTATION  OF  B2  WEN  B2  REPRESENT 

TFE  CATFDDIC  TAFEL  SLOPE 
D3  :  CHARACTER  STRING  USED  TO  FORM  TITl_E2 
D4  :  SAME  AS  02  BUT  FOR  TFE  AMDIC  TAFEL  SLOPE 
LIST  :  OHARACTBR  STRING  USED  FDR  MESSAGES  AND  PRDFPTS 


FEALM  XDAT (2CXX)), YDAT (2000), XDATl (2CX)0), 

+YDAT1(2000) 

REAL*8  M)ATA(2000,2),M)ATA1(2000,2),X1,X2,NB,XY,XI,Y1 
+  ,a)ATA(20OO,2)  ,(3)ATA1  (2000,2) 

INTEGER*2  I  ,  J,N,A1  ,AZ,C1  ,C2,01  ,NI , J1,S1  ,S2,S3,S4,D2,  J2,NB 
INTEGEF»2  A3,A4,CS,C4 

aHPRACTER*8  NArE,NAME2,  TITLE2*32,YNArE*5,LI5T»50,TlTLE3*50 
CHPRACTER*8  D2,D4,D1*5,D3*6 
Tin_E2  = 

YNAME  =  '  ' 

TFE  SCFEEN  IS  CLEAfED  AM)  TFE  FILE  NATE  IS  ENTEFED.  TFE  USER 
CAN  ALSO  EXIT  TFE  PROGRAM  (H  GET  A  LIST  CF  EXI5TIM3  DATA  FILES. 
THIS  PHOrPT  AFTEARS  IN  FIFPLE  WHICH  IS  TFE  OLOR  I RFD  WEMEVER 
A  CHARACTER  FE3=ONBE  IS  EXPECTED  TO  A  PROMPT.  IF  A  FILE  NATE  IS 
ENTEIED  WHICH  DCES  NOT  EXIST  AN  EJHRDR  FESSAGE  WILL  APPEAR  IN  RED 


CALL  aBMDDE(3) 
CALL  aCLEAR(0,5) 


LISl  =  ■  ■ 
cm_  a>CW(l5,12) 

WITE(6,a05)  LIST 
CALL  GLM]V(25,12) 

LIST  =  ENTER  THE  FILE  NATt  TD  BE  LBED' 

WRITE(6,a05)  LIST 
CALL  GOl]W(20,ll) 

LIST  =  A  BLAM<  TO  EXIT,  OR  "DIR”  FOR  A  DIRECTCRY' 

WR1TE(6,805)  LIST 
CALL  GCMCW(36,10) 

REAO(6,800)  NATE 
CALL  a:M:w(w,9) 

IF  (NAME  -EQ.  '  ')  TfEN 

GOTO  900 
EM)IF 

IF  (NAft  .EQ.  DIR'  )  7VEN 
C:ALL  aCLEAR(0,5) 

CALL  a:«J^(l,23) 

PALEE  ENTER  "DIR  *.  /W"  OR  A  BLAN<  TO  CCNTINLJE 
CALL  ajBCR_(Z3,6,l, 1,79,0,5) 

GOTO  15 
ENDIF 

CALL  DATAIN(CDATA,N,NArE,  lER) 

IF  (lER  .ED.  1)  THEN 
(3010  lO 
ENDIF 

Cyil_  DATAM(CDATA,rt)ATA,N,J) 

CALL  CHEn<(«)ATA,N,J,Sl,S2,S3,S4) 

CALL  GBMI}E(3) 

CALL  GnjEAR(0,2) 

J1  =  J 
N  =  N-1 
□1  =  N 
N1  =  N 

TFE  MAIN  rENU  IS  DI31_AYED  IN  QFEEN.  TVE  COLOR  GFEEN  IS  LEED 
^«^E^EV€R  AN  INTEGER  IS  EXPECTED  AS  THE  FESETEE  TD  A  HR>PT. 

IF  TIE  IIWT  IS  NOT  AN  INTEGER  CF  A  CHARACTER  AN  EFUJR  MESSAGE 
IS  DISPLAYED  IN  FED. 

CALL  GI>tX'(15,22) 

LIST  =  '  • 

U«nE(6,805)  LIST 
CALL  00107(15,10) 

LIST  =  1.  FLOr  ORIGINAL  HOTENTIODYNAMIC  CLR^^ 

li«ITE(6,a05)  LIST 
CALL  0007(15,17) 

LIST  =  2.  (3EICRATE  TAFEL  SLOPES  LB1^G  CUBIC  SFLIFE  TEDCI) 

k«ITE(6,a05)  LIST 
CALL  0007(15,16) 

LIST  =  3.  (SSERATE  TAFEL  51_CFtS  UBir«  CENTHA1_  DIFFETENCE  rETTOD 


v 'A'*'.'** 


D-M»  335  DETERHIWITION  OF  TAFEL  CONSTANTS  IN  NONLINEAA 

POLARIZATION  CURVESCU)  NAVAL  POSTGRADUATE  SCHOOL 
NONTEREV  CA  T  E  O'LOUOHLIN  DEC  87 


UNCLASSIFIED 


F/Q  il/'fi.  1 


n  n  n  n 


V«ITE(6,805)  LIST 
CALL  acWW(15,15) 

LIST  =  *4.  USE  flfCnVER  DATA  FILE' 

k«ITE(6,a05)  LIST 

LIST  =  '5.  QMEFLAY  TVD  POTENTICDYNAMIC  CURVES' 

CALL  GCm/(15,14) 
loRITEIAjaDS)  LIST 
CALL  acm/(  15,13) 

LIST  =  '6.  EXIT' 

l•RITE(6,805)  LIST 
CALL  GCm/(29,12) 

LIST  =  'ENTER  1,2, 3,4, 5  OR  6' 
hRITE(6,e05)  LIST 
CALL  GCrCW(39,10) 

FEA0(6,810,EPR  =  22,ia5TAT  =  J4)  M 
CALL  QCMDW(39,9) 

22  IF  ((J4  .hC.  O)  .OR.  (M  .GTT.  6))  TTCN 
CALL  a9fXE(3) 

CALL  GCLEAR(0,4) 

LIST  =  '  ' 

CALL  aCMV(20,17) 

WITE(6,805)  LIST 

LIST  =  'YOU  DID  MJT  ENTER  AN  INTEGER  CF' 

CALL  GCMX'(20,16) 

URITE(6,a05)  LIST 
LIST  =  '\«LLJE  1, 2,3,4  OR  5.' 

CALL  GCM»;(20,15) 

WRlTE(6,a05)  LIST 
LIST  =  'FFESB  ANY  KEY  TO  ONTIMlE:' 

CALL  0(70/(20,13) 
hRITE(6,a05)  LIST 
CALL  aCTO/(39,12) 

(TILL  QIN<EY(I3,14) 

GOTO  20 
EM5IF 

OPTION  1  HAG  BEEN  SLGCTED.  THIS  GE^ERATES  A  GRAPH  AM)  OPTIONAL 
hjdt  of  tve  original  POTENTIODYNAMIC  CLFVE. 

IF  (M  .ED.  1)  IFEN 
DO  25  I  =  1,01 

XDAT(I)  =  M)ATA(I,2) 

YDAT(I)  =  W)ATA(I,1) 

25  (XNTIMJE 

Ml  =  N 
NB  =  O 

YNATE  =  'V/DLTS' 

TnLE2  =  '  ' 

TITLE3='  ' 

CALL  GHARTl(XDAT,YDAT,01,Jl,NArE,XDATl,YDATl,NB,YNArE,TITL£2 
+  ,A1,A2,C1,C2,  TITLES) 
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QOTD  20 
EM)IF 

CFTION  2  HftB  BEEN  SELECTED.  THIS  GE^ERA^ES  TTt  DERIVATIVES  CF 
THE  POTENTIGDVNPrilC  ajR\€;  USING  D-E  CUBIC  SPUirC  lETHDD. 

IF  (M  .EQ.  2)  TVEN 
A1  =  O 
A2  =  O 
Cl  =  O 
C2  =  O 
0  J1  =  J 

N1  =  N 

CflLL  SLDPEl(rt>ATA,Nl,Jl,NPre;,Al,A2,Cl,C2) 

STILL  WITHIN  CFTION  2.  TEE  TTOERflM  OECKS  TO  ENELFE  THAT  TEE 
TAFEL  FEGIGNS  FOR  BOTH  TEE  ANODIC  AM)  CATECDIC  FEGIONB  HA^E 
EEEN  ^1  i-< : 

IF  ((Cl  .GE.  C2)  .CF.  (A1  .QE.  A2))  TEEN 
C/ILL  □Gr>CCE(3) 

CALL  aCLEAR(0,4) 

LIST  =  '  ' 

CALL  acro/(i5,20) 

IF1TE(6,0O5)  LIST 

LIST  =  YOU  DID  NOT  SELECT  BOTH  AN  AMFIC' 

CALL  010/(15,19) 

WRITE(6,905)  LIST 

LIST  *  AM)  CATECDIC  HRPMF  OR  ELSE  TEE!  POINTS' 

CALL  GCrO/(lS,18) 

MRITE(6,805)  LIST 

LIST  =  'EEFE  NOT  PROPERLY  TRAPPED.  IN  ORDER  TO" 

CALL  000(15,17) 

MRITE(6,805)  LIST 

LIST  =  'PREVENT  TEE  PROGRAM  FROM  CRA&IINB  YOU' 

CfiLL  000(15,16) 

WRITE(6,805)  LIST 

LIST  =  'MJBT  START  AGAIN.' 

CALL  CDO(15,15) 
kRITE(6,805)  LIST 
LIST  =  PRESS  ANY  KEY  TO  OFTINLE' 

CALL  000(15,13) 

»i«ITE(6,0O5)  LIST 
CALL  0IM<EY(I3,I4) 

GOTTO  20 
EMOIF 

ALL  fEDUIRENTS  E«ME  BEEN  STTISFIED.  L1^EAR  fEEFESBION  IS 
PERFORED  ON  TEE  TAFEL  FEEIONS  AM)  TEE  FESLL.T5  ARE  GRAFTED 

D1  =  'Bc=  ' 
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VNPTE  =  VOLTS’ 

DO  35  I  =  1,01 

XDAT(I)  =  M3ftTA(I,2) 

YDAT(I)  =  «)ATA(I,1) 

35  GCNTINLE 

X2  =  O 
XI  =  O 
XY  =  O 

NI  =  ABS(C2-C1)  +1 
YI  =  O 
C3  = 

C4  =  J-(44C2) 

DO  40  I  =  C4,C3 

XI  =  XI  r«ATA(I,2) 

X2  =  X2  +  («3ATA(I,2)«*2) 

XY  =  XY  +  (r«ATA(I,2)*M)ATA(I,l)) 

YI  =  YI+  raiATAd.l) 

40  CIKriNLJE 

^e  =  FIDAKNI ) 

81  =  ((XY*XI)-(YI*X2))/((XI*XI)-(NB*X2)) 

82=  ((XY*t«)-(YI*XI))/((X2«0-(XI*XI)) 
YI  =  MDATA(J,1) 

XI  =  ((Y1-B1)/B2) 

X2  =  MDATA(C4,2) 

DX  =  (X2-Xl)/(50.) 

DO  45  I  =  1,50 

YDATl(I)  =  (B2*X2)+B1 
XDATl(I)  =  X2 
X2  =  X2-  DX 
45  OKTIf^ 

82  =  -82 

CPU-  Ca«;(D2,B2) 

X2  =  O 

XI  =  O 

XY  =  O 

^G  =  O 

YI  =  O 

A3  =  J+4+A1 

A4  =  J+4+A2 

DO  50  I  =  A3,A4 

XI  =  XI  +  MDATA(I,2) 

X2  =  X2  +  (f«ATA(I,2)»*2) 

XY  =  XY  +  (rOATA(I,2)m)ATA(I,l)) 

YI  =  YI+  rt)ATA(I,l) 

50  CONTINUE 

NI  =  PB5(A3-A4)+1 
NB  =  FlJOAT(NI) 

81  =  ((XY*XI)-(YI*X2))/((XI*XI)-(NB*X2)) 

82  =  ((XY*NB)-(YI*XI))/((X2«NB)-(XI«XI)) 
YI  =  MDATA(J,1) 


XI  =  ((Yl-Bl)/a2) 

X2  =  t«f)ftTA(A4,2) 

DX  =  (X2-Xl)/(50. ) 

DO  55  I  =  51,100 

YDATKI)  =  (B2»X1)+B1 
XDATl(I)  =  XI 
XI  =  X1+  DX 
CCNTINLE 
D3  =  '  Ba=  ' 

CALL  Cm/(D4,B2) 

X2  =  (XDftTl(500)4-XDATl(501))/2 

TniJE2=Dl//D2//D3//D4 

NB  =  100 

CALL  GRflmi(XDAT,YDAT,01,Jl,Nf*t,XDATl,YDATl,NB,VN««e,T17LE2 
,Al,A2,Cl,C2,TnLE3) 

THIS  PART  flU-OhE  TVE  USER  TO  SELECT  EITVCR  flNDTTCR  ANODIC  CR 
CATVCDIC  TAFEL  REGION  CR  TO  FCVIEW  TTE  GRART  UITHOUT  HAVING 
TO  UETEFWIfE  BOTH  TVE  TAFEL  REGIOE  AGAIN. 

CALL  aCLEAR(0,2) 

LIST  =  '  ' 

CALL  000/(20,16) 
liiRITE(6,a05)  LIST 

LIST  =  *  IF  DE  GRAPH  DID  NOT  SHOW  A  GOOD  CORRELATION' 

CALL  000/(20,15) 

MRITE(6,a05)  LIST 

LIST  =  'YOU  CAN  REPEAT  THE  CALCULATIONS  FOR  THE' 

CALL  000/(20,14) 

MI1TE(6,805)  LIST 

LIST  =  'BRANCH  IN  0LE5TI0N' 

CALL  000/(20,13) 

WRITE(6,805)  LIST 

LIST  =  '1.  REPEAT  CALCULATIONS' 

CALL  000/(20,11) 

URITE(6,805)  LIST 

LIST  =  '2.  EXIT  THIS  PORTION' 

CALL  000/(20,9) 

hRITE(6,805)  LIST 

LIST  =  ENTER  A  1  OR  A  2' 

CALL  000/(20,7) 

WRITE(6,0O5)  LIST 
CALL  CDO/(25,6) 

READ(6,B10,EPR=56)  M 
IF  (M  .ED.  1)  TPCN 
GOTO  30 
EM)IF 

IF  (M  .EQ.  2)  TIEN 


EMJIF 

EM}IF 

T>E  CENTRRL  DIFFEFNCE  fCTHD  IS  TO  BE  USED.  Tl-E  OFTirN  FfOEEDS 
IN  TFE  SPt£  MPMSER  flB  CPTIGN  2 

IF  (M  .EQ.  3)  IVEN 
A1  =  O 
A2  =  O 
Cl  =  O 
C2  =  O 

CWLL  SLOFEIMMTA.QI , J1,NP^E, A1 , A2,C1  ,C2,S1  ,S2,S3,S4 ) 

IF  ((Cl  .GE.  C2)  .CR.  (A1  .(BE.  A2) )  TVEN 
CPLL  asro£(3) 

CALL  aCLEAR(0,4) 

LIST  =  '  ' 

(>«_L  CD>O/(lS,20) 

WRITE(6,805)  LIST 

LIST  =  'YOU  DID  NOT  SELECT  BOTH  AN  ANODIC' 

CALL  (rr<CV(15,19) 

|pRITE(6,a05)  LIST 

LIST  =  'AM)  CATHODIC  BRANCH  OR  ELSE  TVE  POINTS' 

CALL  (Z7CV(15,18) 

WRITE(6,805)  LIST 

LIST  =  'dCIE  NOT  PFCPEFL.Y  TRAPPED.  IN  ORDER  TO' 

CALL  GDCWI  15,17) 

UFITE(6,a05)  LIST 

LIST  =  'PFE^ENT  TTC  PROGRAM  FF01  CRABHIM3  YOU' 

CAUL  aCrfX'(15,16) 

WRITE(6,a05)  LIST 

LIST  ==  'MUST  START  A&UN.' 

CALL  0110/(15,15) 
tiiRITE(6,805)  LIST 
LIST  =  'FPE55  ANY  KEY  TO  OCNTIMJE' 

CALL  (X10;(15,13) 

URITE(6,a05)  LIST 
CALL  aiN<EY(I3,I4) 

GOTO  20 
EM)IF 
Dl='Bt=  ' 

YNAME  =  'VOLTS' 

DO  60  I  =  1,01 

XDAT(I)  =  M)ATA(I,2) 

YDAT(I)  =  M)ATA(I,1) 

CONTINUE 
X2  =  O 
XI  =  O 
XY  =  O 

NI  =  ABB(C2-C1)  +1 
YI  =  O 

C3  =  J-(Z+C1) 


C4  =  J-(24C2) 

DO  7X>  I  =  C4,C3 

XI  =  XI  +  rOATA(I,2) 

X2  =  X2  +  (MQATA(r,2)»*2) 

XY  =  XY  (r«ATA(I,2)««ATA(I,l)) 

YI  =  YI+  rOATA(I,l) 

CDNTIMJE 
NB  =  FlJQAT(NI) 

B1  =  ((XY*XI)-(Yr*X2))/((XI*XI)-(Ne*X2)) 
H2  =  ((XY*NB)-(YI«XI))/((X2»e)-(XI*XI)) 
Yl  =  K)ATA(J,L) 

XI  =  ((Y1-B1)/B2) 

X2  =  ««ATA(C4,2) 

DX  =  (X2-Xl)/(50.) 

DO  80  I  =  1,50 

YDATl(I)  =  (e2*X2)+fll 
XDATKU  =  X2 
X2  =  X2-  DX 
GCNTINLE 
H2  =  -82 
CALL  CO«/(D2,B2) 

D3  ='  Ba=  • 

X2  =  O 

XI  =  O 

XY  =  O 

NS  =  O 

YI  =  O 

A3  =  J+2^1 

A4  =  J+2+A2 

DO  90  I  =  A3,A4 

XI  =  XI  +  «)ATA(I,2) 

X2  =  X2  +  (ra}ATA(I,2)««2) 

XY  =  XY  +  (rCATA(I,2)*M]ATA(I,l)) 

YI  =  YI+  MJATA(I,1) 

CONTINUE 

NI  =  AB5(A2-A1)+1 
^e  =  RJOAT(NI) 

81  =  ((XY*XI)-(YI»X2))/((XItXI)-(r6»X2)) 

82  =  ((XY««)-(YI»XI))/((X2»e)-(XI*XI)) 
Yl  =  rCIATA(J,l) 

XI  =  ((Y1-B1)/B2) 

X2  -=  rC)ATA(A4,2) 

DX  =  (X2-Xl)/(50. ) 

DO  lOO  I  =  51,100 

YDATl(I)  =  (B2*X1)+B1 
XDATl(I)  =  XI 
XI  =  X1+  DX 
COYTIMjE 
CALL  CaN^(D4,B2) 

TITLE2=01//D2//D3//D4 
NB  =  100 
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CflLL  a=WPHl(XDAT,ynAT,Dl,Jl,NPfC,XDATl,YDATl,NB,YMfirC,TITUE2 
,Al,A2,ni,C2,TlTLE3) 
mLL  UCLJ£PR(0,2) 

LIST  =  •  ' 
cnu-  Gcm/(20,i6) 

WRITE(6,805)  LIST 

LIST  =  '  IF  TTE  SWH  DID  ^fJT  SHEW  A  GOD  CCFRELATIGN' 

CflLL  aC«W(2D,15) 
li«nE(6,a05)  LIST 

LIST  =  'YOU  CAN  REPEAT  TFE  CAUCLLATIONB  FCR  TVE' 

CALL  CCrClV(20,14) 

WRITE(6,a05)  LIST 

LIST  =  ‘BRPNm  IN  OLESTION' 

CALL  GCrCW(20,13) 

WRI7E(A,a05)  LIST 

LIST  =  '1.  FEPEAT  CALCLLATIONS' 

CALL  acrtW(20,ll) 

W9I1E(6,80S)  LIST 

LIST  =  '2.  EXIT  THIS  PCRTiaN' 

CALL  acrfV(20,9) 

URITE(6,805)  LIST 
LIST  =  'EMIER  A  1  OR  A  2' 

CALL  □□-CW(20,7) 

WRlTE(6,a05)  LIST 
CALL  acrtW(25,6) 

READ(A,810,ERR=s5e)  M 
IF  (M  .ED.  1)  TPEN 
GOra  57 
6M)IF 

IF  («  .ED.  2)  THEN 
GOTO  20 
ELSE 

GOTO  58 
ENDIF 
EMJIF 

CFTICfl  4  RETXJFN5  TEE  USER  TO  TFE  BEBIWINB  OF  TVE  PRDGRAfI 

IF  (M  .EQ.  4)  TVEN 
GOTO  lO 
EM)IF 

CFTICN  5  IS  TVE  CXEFEAY.  CNlY  TVE  PaTENTIODYNPniC  CLF\ES  CAN 
BE  STAfVED  OR  PLOTTED-  TVE  OTIGINAL  FILE  IS  STILL  AMAIL/^LE 
FOR  FURTHER  MANIPLLATICW 

IF  (M  -EQ.  5)  TVEN 
TriLE3='  • 

CALL  GCLEPR(0,5) 

LIST  =  '  ' 

CALL  QCm/(  15,20) 
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WRI7E(&,a05)  LIST 
CfLL.  aCMCW(25,12) 

LIEiT  =  ENTER  TEE  FILE  NPfE  TD  BE  USED' 

V«I7E(6,a05)  LIST 
CflLL  aCMCN(20,U) 

LIST  =  A  a-PM<  TO  EXIT,  CR  "DIR"  FTP  A  DIRECTORY' 
i«PnE(6,a05)  LIST 
GflLL  001^(36,10) 

READ(6,aOO)  NPrE2 
CALL  aCMDW(^,9) 

IF  (NArE2  .ED.  '  '  )  TTEN 

GOTO  21 
END  IF 

IF  (NPfE2  .ED.  DIR  )  TTEN 
OEL  aCLEAR(0,5) 

LALC  aC«V(l,23) 

PAUSE  ENTER  “DIR  *.  /W“  OR  A  BLAMC  TD  OONTINE' 

CALL  aLJSCFa-.(23,6, 1,1,79,0,5) 

□OTD  105 
EMTIF 

DO  750  I  =  1,01 

XDAT(I)  =  «>ATA(I,2) 

YDAT(I)  =  MDATA(I,1) 

CETTriME 

CALL  DATAIN(CDATAl,a2,NArE2,IER) 

IF  (lER  .ED.  1)  TTEN 
GOTO  105 
EM)  IF 

CALL  DATAn(CDATAl,MDATAl,a2,J2) 
cm.  OECK(MDATAl,a2,J2,Sl,S2,S3,S4) 

DO  760  1=1,02 

XDATKI)  =  M)ATA1(I,2) 

YDATl(I)  =  M)ATA1(I,1) 

CENTIME 

TITTJE2  ='C]^EFLAY  '//NATE//*  ' //NPrE2 
YNATE  =  'VOLTS' 

NAfE2=NArE 
NATE  =  '  ' 

CALL  SWHl(XDAT,YDAT,01,Jl,NPrE,XDATl,YDATl ,02, YNATE, TITLE2, 
A1,A2,C1,C2,TITI£3) 

NATE  =  NPfE2 
CALL  aCLEAR(0,2) 

GOTO  21 

ENDIF 

IF  (M  .ED.  6)  TTEN 
GOTO  900 

nsF 

GOTO  20 
EM)IF 
Ftf»1AT(AB) 

FOR1AT(A50,\) 
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THIS  aJBRDUTIlvC;  REALS  TVE  DATA  FILE  GEJEFTTATED  BY  fl-E 
nXEL  351  UHICH  HAS  SEN  THANEFEFEED  LNTO  A  DISK 


□DATA 


=  TVE  ARRAY  WHICH  CXNTAININB  POTENTIAL  AND 
OJFfENT  SNSI7Y 
=  TVE  NATE  S  TVE  DATA  FILE 
=  CHARACTER  STRINB  FTESENT  IN  TVE  DATA  FILE 
=  CHARACTER  STRING  PRESENT  IN  TVE  DATA  FILE 
=  NLIEER  CF  DATA  POINTS  IN  TVE  FILE, 
DETEFWIfED  DURIM3  TVE  PROGRAM 


SLJGROLrTIS  OATAIN  (CDATA,N,NAVE,  lER) 

FEALta  CDATA(2CXX>,2) 

INTEaER«2,  N,I 

CHARACTER  Cl*17,CZM),NArE*a,LIST*50,C3*50 

□PEN  (INIT  =  10,  EFR  =  200,F1LE  =  NAME,  STATLE  =  CLD  ) 

DO  50  N=l,2000 

fEAD  (10,900,ET«=60)  Cl,CnATA(N,l ) ,C2,CDATA(N,2) 
CCNTINLE 
SMINO(IO) 

DO  70  I  =  1,2000 

REA0(10,«,EFR  =  100) 

CCNTIME 

IF  (I  .NE.  N)  TVEN 
CALL  QSM0DE(3) 

CALL  aCLEAR(0,4) 

LIST  =  '  ' 

CALL  aCNCW(15,20) 
l«RITE(6,905)  LIST 

LIST  =  AN  EFROR  HAS  BEEN  ENCOLNTEFED  IN' 

CALL  aCM0V(15,19) 

WRITE(6,905)  LIST 

LIST  =  TIE  FCFriAT  CF  TVE  DATA  FILE.' 

CALL  aa*tV<15,18) 

WRITE(6,905)  LIST 

LIST  =  THIS  USUALLY  IFNCLVES  A  SWITCH  BETWEEN' 

CALL  acro/(i5,i7) 

WRITE(6,905)  LIST 

LIST  =  FLOATING  POINT  AW  EXPOENTIAL  NJTATICN' 

CALL  aCMCN(15,16) 

WRITE(6,905)  LIST 
LIST  =  AROUM)  ECOCR.' 

CTEL  QCm/1 15,15) 

UR1TE(6,905)  LIST 

LIST  =  TO  COFFECT  THIS  FRCBLEM  EXIT  TVE  PROGRAM' 
CALL  0CM0W(15,13) 

WFITE(A,905)  LIST 

LIST  =  AM)  UEE  TVE  DOS  "EDLIN”  FEATUE.  TVE' 


CALL  GDtWC  15,12) 

WRITC(6,905)  LIST 

LIST  =  'PBCBLBI  BEGINS  AT  LINE  ' 

CAUL  tX3«W(15,ll) 

WRI  fE(6,905)  LIST 
CALL  aC10/(3B,ll) 
l«RITE(6,915)  N 

LIST  =  PFESS  ANY  KEY  TO  CCMriMJE' 

CALL  COCV(lS,8) 
l««ITE(6,905)  LIST 
CALL  aiM<EY(I3,I4) 

GOTO  990 
EM)IF 
GCrra  990 
CALL  a5rO)E(3) 

CALL  aCLEAR(0,4) 

LIST  =  •  ' 

CALL  OC — ■  (15,20) 

WRnE(6,905)  LIST 

LIST  =  Tl-E  FILE  NATE  EMTEFED  DOES  (OT  EXIST.' 

CALL  aCMaV(15,19) 

WRnE(6,905)  LIST 

LIST  =  TO  EXAm^E  TVE  DIRECTORY  (T  POSSIBLE  FILES' 
CALL  (30*0/(15, 18) 
hRI7E(6,905)  LIST 

LIST  =  ENTER  "DIR"  MEN  ASKED  TO  ENTER  TTE  FILE' 
(^U_  (30*0/(15,17) 

MTITE(6,905)  LIST 

LIST  =  I^IArE:  RXLflED  BY  THE  DOS  CTTMAND. ' 

CALL  acro/(&,905) 

WRnE(6,905)  LIST 
LIST  =  '"DIR  *.  /W"' 

CALL  (30*0/(35,16) 

UFnE(6,905)  LIST 

LIST  =  'FTESS  ANY  KEY  TO  CINTIISUE' 

CALL  acm/(15,14) 

WRITE  (6,905)  LIST 
lER  =  1 

CALL  aiM<EY(I3,I4) 

(SOTO  990 

F0R1AT  (A17,GB.4,A7,G10.4) 

FCR1AT(A50,\) 

FORMAT ( I 4, \) 

(ETLFN 


THIS  SLBRUTir^  MPNIPUL^TES  THE  DATA  IN  TVE  ARRAY  CDATA 
TVE  PREFER  FORM  AW)  IDTENIFIES  ETTFR 

CDATA  =  TVE  ARRAY  CTKTAININB  POTENTAIL  AW)  CURRENT 
DENSITY 

WJATA  =  CDATA  CIFfERTED  TO  POTENTIAL  AND  UDG 
CURRENT  I»GITY 

N  =  NUMBER  OF  DATA  POINTS  IN  TVE  DATA  FILE 

J  =  LDCATIGN  WITHIN  rt)ATA  WHERE  EECFR  IS  LOCATED 

I  =  CO  LOCP  CCLNTER 

XI  =  ^fflSOUJTE  V/ALIE  CF  CUFFENT  DENSITY  AT  I 
USED  FOR  CETER1ININB  ECLRR 
X2  =  ABSCLUTE  VALUE  OF  CURRENT  DENSITY  AT  I-l 
USED  FOR  DETER1ININB  ECORR 
X3  =  UJGIO(XI) 

SUEROUrirE  DATAM  (CDATA,W)ATA,N,J) 

REAL*a  rt)ATA(20OO,2)  ,CDATA(2000,2) 

INTEGER«2  N,J 
X2  =  100000. 

OOWERT  CLRFENT  DENSITY  INTO  LOG  CURENT  DENSITY 

DO  100  I=l,N-2 

XI  =  ABS(CDATA(I,2)) 

X3  =  LOGIOIXI) 

IF/TVEN  TEST  USED  IN  DETEFMINIW3  LOCATICN  CF  ECDDR 

IF  (  XI  .LT.  X2)  TVEN 
J  =  I 
X2  =  XI 
ELSE 

X2  =  X2 
EW)IF 

MDATA(I,1)  =  CDATA(I,1) 

MDATA(I,2)  =  X3 
3  CCNTIWE 
FETUPR 
EW) 


a«>  m-H  «—  »»A  jn  ^*»<i  •■**•»**• 
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7HIS  aUEFajTltE.  PEFFCmS  7WD  FUsCTIUe.  ITS  FIFBT  FUNCTICW 
IS  TO  O-ECK  TO  ENBLFE  THAT  TFE  POTENTIAL  VmJES  RECORDED 
ACCURATELY  FCFIJECT  TVC  ACTUAL  POLARIZATION.  Tl-E  SECCN}  IS 
TO  DETER1I^E  TVC  NLTCER  OF  POINTS  IN  TVE  AMJDIC  AND  CArHDDIC 
BRANCHES  WHICH  MAY  IE  ttm  IN  THE  CENTRAL  DIFFEPNCE  METHOD. 

M)ATA  :  THAT  ARRAY  CCNTAINIIVC  POTENTIAL  AM)  LOB  CLPRENT 
DENSITY 

N3  ;  THE  NUMBER  OF  POINTS  IN  TVE  CATVHODIC  BRANCH  UHICH 
CAN  BE  IFEn  IN  TVE  CENTTV^  DIFFEFENCE  rETHOD 
N4  :  TVE  fOEER  OF  POINTS  IN  TVE  ANODIC  BRANCH  l«HICH 
CAN  BE  I  EFT)  IN  TVE  CENTFJAL  DIFFEREE  rETVO) 

OTVER  MARIBLES  ARE  tm)  tERELY  IN  OECKING  REGUIFED  CONDITICIE 

SUmOUTirE  CVECK(MDATA,N,J,N1,N2,N3,N4) 

FEAL«8  ra)ATA(2000,2) 

INTEGERS  I,J,N,N1,N2,N3,N4 
CHARACTER*50  LIST 

TVE  FIRST  INOEMENT  IN  POTENTIAL  FOR  TVE  CATVCDIC  BRAMH  IS 
DETERMINED 

Y  =  rtDATA(J,l) 

Y1  =  M)ATA(J-1,1) 

Z1  =  ABB(ABS<Y1)-^«S(Y)) 

Z1  =  ANINT(Zl/O.CX»l) 

SUBSOUENT  POTENTIAL  INCFEMENTS,  VRICH  MJBT  BE  EQUAL  TO  DE 
FIRST  ARE  OEEKED.  LOG  CLHENT  DENSITY  VALLES  ATE  ALSO  OECKED 
TO  ENSURE  THAT  TVEY  WILL  MIT  fESLLT  IN  TVE  CWLCULATICN  OF  A 
SLOPE  OF  IVFINITE  VALLE. 

DO  20  I  = 

Y  =  M)ATA(I,1) 

Y1  =  rf)ATA(  1-1,1) 

X  =  M)ATA(I,2) 

XI  =  M)ATA( 1-1,2) 

Z  =  ABS(ABS(Y1)-ABS(Y)) 

Z  =  ANINT(Z/0.0001) 

GJE  OF  TVE  TWO  FEEDIFED  CONDITIONS  HA6  FAILED. 

IF  ((Z  .^E.  Zl)  .CR.  (  XI  .LE.  X  ))  GOTO  25 

CDNTIVUE 
Ml  =  1+2 


& 

i 


THE  FIRST  INCFEl-ENT  IN  POItNTIflL  FDR  TVE  P^a5IC  BRPM>J  IS 
DETERMI^ED 

X  =  rt)ATA(J+l,l) 

Y  =  M}ATA(J+2,1) 

Z1  =  #«B(X-Y) 

11  =  ANIMTIZl/O.OOOl) 

SI HBFIUFT^fr  RJTENTIAL  IMTO^NrS,  l«HICH  rUCT  BE  EQLiPL  TO  TIE 
FIRST  flfE  OECKED.  I_CE  CURRENT  IHEITY  VCLLES  AFE  ALSO  OECKED 
TO  ENSURE  THAT  TVEY  WILL  NOT  FESULT  IN  TTE  CALCULATION  CF  A 
SLIFE  OF  I^FINI7E  VALJUE. 

DO  30  I  =  J+2,N-1 

Y  =  rt)ATA(I,l) 

Y1  =  MJATAC 1+1,1) 

X  =  rt)ATA(I,2) 

XI  =  MJATAI 1+1,2) 

Z  =  ABB(Yl-Y) 

Z  =  ANINTIZ/O.OOOl) 

GTE  OF  TVE  TWO  REOUIFED  COMIITIONB  NAB  FAILED. 

IF  ((Z  .IE.  Zl)  .OR.  (  XI  .LE.  X  ))  GOTO  35 

CONTIME 
N2  =  I  -  2 
N3  =  J  -  N1 
N4  =  N2  -  J 

TVE  rOEER  OF  POINTS  AMAlLAOE  FDR  EAOT  ERAfCH  IS  COtWED 
WITH  A  MIMMILfl  NJEER  CF  15.  IF  lEITVER  BRMEH  HAS  ENOUGH 
AMAILABLE  POINTS  VARIOJB  lEBSACES  ARE  DI9EAYED  LISTIIE 
PCBSIBLE  CAUSES  AM)  CPTIGN5. 

IF  ((N3  .LE.  15)  .OR.  (W  .LE.  15))  TVEN 
CALL  □CLEAR(0,4) 

LIST  =  •  • 

CALL  (I>CN(20,20) 

WRITE<6,805)  LIST 

LIST  =  TVE  CENTTVE  DIFFEREE  lETVCD  FEQUIRES  THAT' 

CALL  GC«W(20,19) 

WRITE(6,a05)  LIST 

LIST  =  TIE  dFVENT  EE  lEASJED  AT  EQUAL  DIFFERENTIALS- 
CALL  aDEV(20,18) 

MTlTE(6,a05)  LIST 

LIST  =  OF  TIE  AREIED  \4a_TAGE.  TVE  PAR  rtBEL  351- 
CALL  GCm/(20,17) 
l««ITE(Z,,0O5)  LIST 

LIST  =  RECOFDS  TIE  CLFRENT  AT  IfCFEMENTS  CF  2.0  aH' 

CALL  000/(20,16) 
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»««ITE(6,a05)  LIST 

LIST  =  ’UNTIL  TT-E  CyVt43E  IN  tURFENT  EXCEEDS  A  PPEBCJIIBED' 
cflLL  acro/(20,i5) 
hRrTE(6,a0b)  LIST 

LIST  =  LEVEL.  AT  THAT  POINT  TVE  351  BEGINS  RECXFDINB’ 

CPU-  acn:v(20,i4) 

»««nE(6,a05)  LIST 

LIST  =  '7TE  OLFFENT  AT  ItEFEJENTS  OF  0.5  irW.  AT  PN’ 

CPU-  CD1W(20,13) 

WRnE(6,a05)  LIST 

LIST  =  HIGKER  PFELIED  VO-TPGE  TVE  OHPNEE  IN  CURRENT' 
CPU-  aC«W(20,12) 

»«RITE<&,a05)  LIST 

LIST  =  'DECREASES  TO  A  POINT  WVEIE:  TVE  CURRENT  15  AGAIN' 
CPU.  QCraV'(20,ll) 

MRlTE(6,a05)  LIST 

LIST  =  'fEASURED  AT  IfCREMENTS  CF  TVE  APR-IED  VOLTAGE' 
CALL  acriv(2o,io) 

»«RITE(6,805)  LIST 
LIST  =  'EQUAL  TO  2.0  nN' 

CALL  acrtw(2o,9) 

UfTITE(6,a05)  LIST 

LIST  =  PRESS  ANV  KEY  TO  CCNTINLE' 

CPU-  acrcN(20,7) 

MTITE(6,a05)  LIST 
CALL  aiM<EY(I,MB) 

CALL  aCLEAR(0,4) 

LIST  =  '  ' 

CALL  acrO/(20,20) 

URI1E(6,805)  LIST 

LIST  =  THE  DATA  FECDRDED  FOR  THIS  EXPERIiENT  EITHER:' 
CALL  acrtN(20,19) 

WRITE(6,805)  LIST 
LIST  =  '  ' 

CALL  aCWV(6,18) 

WTITE(6,805)  LIST 

LIST  =  1.  REFLECTS  THAT  TVE  CHAMFB  IN  CUFFENT’ 

CALL  acrOM(20,17) 

V«ITE(6,805)  LIST 

LIST  =  'rEASJEIENT  DESCRIBED  EARLIER  DID  NOT  OCCLR' 

CALL  0CTtN(20,16) 

UfTITE(6,0O5)  LIST 

LIST  =  AT  A  POINT(S),  RELATIVE  TO  ECOCF,  THAT  WOLLD' 

CALL  GCrCN(20,15) 

WRITE(6,805)  LIST 

LIST  =  GENERATE  ENOUGH  DATA  POINTS  FOR  THE  UGE  OF  TVE' 
mu.  QCWV(20,14) 

MTITE(6,a05)  LIST 

LIST  =  CENTRAL  DIFFERENCE  REIHOD. ' 

CALL  000/(20,13) 

MTITE(6,805)  LIST 


I 
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CPLl-  aCICW(20,7) 

WRITE(6,a05)  LIST 

LIST  =  FWESS  #?NY  KEY  TD  COMTIMJE' 

CflLL  aCrt3V(20,5) 

WRnE(6,S05)  LIST 
CflLL  QIM<EYn,MB) 

Cm_  aCLEflR(0,2) 

LIST  =  1.  DO  rOT  OJWtLT  TT-E  PFBBLEM. ' 

CflLL  acm/(25,13) 

MRITE(6,a05)  LIST 

LIST  =  '2.  TRY  ««  CORRECT  TVE  PROBLEM.  ' 

CflLL  010/(25,11) 

WRITE(6,a05)  LIST 

LIST  =  'ENTER  A  1  OR  A  2' 

CPLL  CDO/(32,9) 
fcFITE(6,a05)  LIST 
CALL  GICm/(39,8) 

)EAD(^>,810,ER1  450,  lOSTAT  =  J4)  M 

IF  ((M  .GT.  2)  .OR.  (J4  .ME.  O) )  THEN 
CALL  Oa_EAR(0,4) 

LIST  =  •  • 

CALL  000/(21, 16) 
l«I7E(6,805)  LIST 

LIST  =  YCU  DID  NOT  ENTER  AN  INTEGER  CF' 

CALL  a>0/(20,15) 

WRITE(6,a05)  LIST 
LIST  =  'VALLE  1  CR  2' 

CALL  030/(37,14) 

WRI7E(6,805)  LIST 

LIST  =  'PRESS  ANY  KEY  TO  CONTINUE' 

CALL  000/(25,12) 

WRITE(6,805)  LIST 
CALL  (300/(39,11) 

CALL  OIM<EY(I,MB) 

GOTO  420 
EM)IF 

IF  (M  .ED.  1)  THEN 
N?  =  1 
GOTO  1000 
EMTIF 

AN  ATTOET  IS  BEING  MADE  TO  IXFWECr  TVE  POTENTAIL  VALLES  SUCH 
THAT  TTEY  ACDLRATELY  REFLECT  TVE  EXPERITENTAL  FDTENTIALS. 
BLOCKS  500  AM)  600  COTFECT  TVE  PROBLEM,  BUT  ON_Y  OVE  OF  TVE 
TWO  BLOCKS  ARE  USED  FCR  ANY  GIVEN  FILE. 

IF  (M  .ED.  2)  TVEN 
DO  500  I  =  1,N,2 
Y1  =  M)ATA(I,1) 

Y2  =  M)ATA(I+1,1) 

IF  (Y1  .ED.  Y2)  TVEN 


W}ATA(I,1)  =  t«I)ATA(I,l)  -  0.000500 
N9  =  1 
EM)IF 

500  CDNTINLE 

DO  600  I  =  2,N,2 
Y1  =  rOATAd,!) 

Y2  =  rCATAC  1+1,1) 

IF  (Y1  .ED.  Y2)  7V€N 
rOATAd,!)  =  ^f>ATAd,l)  -  0.000500 
W  =  1 
EM)IF 

600  cxmir^LE 

IF  (N9  .ED.  1)  TFEN 
CflLL  aCLEPR(0,2) 

LIST  =  '  • 

CflLL  aCMDW(30,14) 

WRITE(6,a05)  LIST 

LIST  =  COHRELTICN  SLCCESSRL- 

Q5LL  acn3/(30,13) 

li«nE(6,a05)  LIST 

LIST  =  LEE  EITVCR  ItTVOD' 

CALL  aCWV(30,ll) 
li«nE(6,0O5)  LIST 
LIST  =  PRESS  ANY  KEY  TD  CONTINUE’ 

CALL  acro/(3o,9) 

(««ITE(6,a05)  LIST 
CALL  QIM<EYd,r«) 

GOTO  10 
ELSE 

CALL  aCLEAR(0,4) 

LIST  «  'CDFreETION  UCHXESSFUL  ’ 

CALL  aCM[V(30,13) 

V«nE(6,805)  LIST 

LIST  =  'LEE  EITICR  CLBIC  SFLIft  METKI) 

CALL  aCMCW(30,ll) 

taRITE(6,0O5)  LIST 

LIST  =  'PRESS  ANY  KEY  TO  CONTII^’ 

CALL  amDV(30,9) 

UFnE(6,a05)  LIST 
CALL  aiM<EYd,MB) 

GOTO  1000 
EMJIF 
EM)IF 
eM)IF 

905  F0FWAT(A5O,\) 

010  FERMAT  (ID 

1000  END 


XDAT  :  T>€  X  VWUJES  CF  TVE  FIRST  GR«=H 

YDAT  :  THE  Y  VALLES  OF  TFE  FIRST  (3WH 

XDATl  :  TVE  X  VRLIES  CF  TIE  SEEOND  GRAW 

YDATl  :  TVE  Y  VALLES  CF  THE  SECCWD  GRAPH 

TX  :  XDAT  AS  ORIGINAL  FWSSED  AM)  SCALED.  UEED  TO  ALUOW 
TVE  USER  TO  RESCALE  TVE  AXIS  AM)  MAINTAIN  ALL  TVE 
VALLES  ORIGINALLY  PASSED 
TV  ;  USED  AS  TX  BUT  FOR  YDAT 

TXl  ;  USED  AS  TX  BUT  FOR  XDATl 

TYl  ;  USED  AS  TX  BLIT  FDR  YDATl 

IN  mis  SUBRCUTirE  ALL  ROPTS  AfE  FEPLIED  TO  BY  TVE  LEE  OF 

SLFT  KEYS.  THIS  KEANS  THAT  A  CARRIAGE  RETIFW  IS  NOT  USED. 

SUBEUTlfE  GRAFHl (XDAT,YDAT,N3,J,NAtE, XDATl, YDATl  ,Nft  , T I TLEZ 

+A1  ,A2,C1  ,C2,  TITLES) 

fEAL*4  XDAT(2000)  , YDAT ( 20CX) ), XDATl (2000)  .YDATKZOOO)  ,F,F1  ,F2 
fEAL*4  TX(2000),TY(2000),TXl(2000),TYl(2(X»),X2(l),Y2(i) 

INTEBER*2  K2,K3,K5,IC6,L,XL,XR,YL,YR,MB,A1  ,A2,C1  .LS.M,  lA(aO) 
INTEGER*2  I , J ,N,K7,K8,K9,  JCOL(IO)  ,JR0W(lO)  ,N1,NB,N3,N4,N5,N6,N7 
IN1EGER42  L1,L2,L3,L4 

CHARACTER  NArE*8,YNAME*5,FACrT*6,XNAKE*12,TnLE*15, TITL£2*32 
CHARACTER*28,  GLES,  TITLE3*50 
N7  =  1 
M  =  O 

N  =  MAX(N3,N4) 

N3  =  N3  -  3 
N4  =  N4  -  3 
IF  (N4  .LE.  O)  N4  =  1 

TVE  AFFtAYS  AFE  SCALED  SUDT  THAT  TVE  Y  AXIS  VALLES  AFE  DISCERNIBLE 
FROM  EACH  OTVER 

F  =  MIN(ABS(YDAT(1)),ABS(YDAT(N3))) 

FI  =  MIN<ABS(YDAT1(1)),#«S(YDAT1(N4))) 

IF  (N4  .GT.  2)  TVEN 
F2  =  MIN(F1,F) 

ELSE 

F2  =  F 
EM)IF 

IF  (F2  .GE.  1.)  TVEN 
F2  =  1.0 
FACT  =  x  1- 
EM)IF 

IF  ((F2  .GE.  .1)  .AM).  (F2  .LT.  1.))  TVEN 
F2  =  10 
FACT  =  X  10' 

EMIIF 
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IF  ((F2  .QE.  .01)  .PM).  (F2  .LT.  .1))  rV€N 
F2  ^  100. 

FPCT  =  X  100’ 

QSCIF 

IF  ((F2  .GE.  .001)  .PM).  (F2  .LT.  .01))  THEN 
F2  =  1000 
FPCT  =  'x  1000' 

ELSE 

IF  (F2  .L£.  0.0010)  T>€N 
F2  =  1000 
FPCT  =  x  1000* 

EM)1F 

ENDIF 

DO  6  I  =  1,2000 

YDAT(l)  =  F2*YDAT(1) 

YDPTKI)  =  F2*YDPT1(I) 

CCNTIMJE 
NO  =  N3 
N6  =  N4 

TVe  PEFMflfCNr  PmPYS  PfE.  FILLED 

DO  10  I  =  1,2000 
TX(I)  =  )(DAT(I) 

TY(1)  =  YDPT(I) 

TXl(l)  =  XDATl(I) 

TYl(I)  =  YDATl(I) 

CINTIME 

mE  GRPFH  AXIS  PND  WIMXIW  PRE  DETERilNED 

XmX  =  -1000 

XMIN  =  1000 

VWX  =  XMAX 

YMIN  =  XMIN 

XMAXl  =  XMAX 

XMINI  =  XMIN 

VMAXl  =  XMAX 

VMINl  =  XMIN 

JCOLl  =  90 

JCDL2  =  610 

JR>»1  =  35 

JFOO  =  300 

DO  200  I  =  1,N3 

YMAX  =  MAX(YDAT(I),YMAX) 

XMAX  =  MAX(XDAT(I),XMAX) 

XMIN  =  MIN(XDAT(I),XMIN) 

YMIN  =  MIN(VnAT(I),VMIN> 

(XNTIMJE 

DO  210  I  =  1,N4 

XMAXl  =  f1AX(XDATl(I),XMAXl) 


VMftXl  =  MAX(YDATi(I),yMAXl) 

XMINl  =  mN(XDATl{l),XMINl) 

VT1IN1  =  MIN(yDAri(  n.YMINl) 

CIWTIf4£ 

YFIN  =  (MAX(YmX,ymXl) )  +  .1 

VST  =  (MIN(yMIN,YmNl))  -  .1 

XST  =  (niN(  XniN, XMINl ) )  -.1 

XFIN  =  (MftX(XMftX,XMftXl))  +  .1 

YMIN  =  VST 

YMftX  =  YFIN 

XMftX  =  XFIN 

XMIN  =  XST 

YCRE  =  YMIN 

XCR3  =  XMIN 

XNPME  =  ijaG  I  A/cm2' 

TITLE  =  •  POTENT  ICDVNPMIC' 

ICPT  =  O 
YC1^«RX  =  0.3 
ASPECT  =  1.5 

TFE  GRPm  AXIS  PRE  DRAMM  AM)  LAffl  FO.  TITLES  ARE  PLACED 
CALL  a91XE(16) 

CA_L  a=LDT(JCCLl,JCnJ2,JFa«»l,JRrW2,XMlN,XMAX,VMIN,YMAX,XaRG,YCRG, 
♦  ICPT ,  VCMERX  .ASPECT ) 

XMAJCR  =  (XFIN-XST)/10. 

XMirOR  =  O 
LABEL  =  -I 
l«EC  =  2 

VmJCR  =(VFIN-YST)/10. 

YMirCF  =  O 

N<NT  =  5 

MXJTS  =  O 

IKOjCR  =  1 

ISVrHl_=-2 

KLRSYn=l 

CALL  0CCLCR(0,0) 

CALL  aBETlF(^iXJT5,IKrlJJT,ISV^B□L,KLRSVM) 

CALL  aXAXIS(XST,XFlN.XI>V\ICR,XMir«CJR,LABEL,M]EC) 

CALL  OYAXIS(YST,VFIN,YMAJCR,VMIfCR,LA«EL,M]EC) 

CALL  QXTICS(tM<NT,jai-,JFO«») 

CALL  aETXT(12,XNPrE,2,340,1.0) 

CALL  OYTICSINKNT.JCDL.JRUW) 

CALL  OSTXT ( 5, YIMATt, 2, 1,246.-1) 

CALL  0GTXT(6,FACT, 2, 1,162,-1) 

CALL  QBTXTI 15, TITLE, 2,254,335,0) 

CALL  QDTXT(B,NPrE,2,3B2, 335,0) 

CALL  nBTXT(32,TITLE2,2,254,321,0) 

CALL  aSTXT(50,TITLE3,2, 182,307,0) 

IKULCR  =  2 

CALL  fEETlJ^(|vaJTS,II<ajTJ,ISYMBa.,KLRSYM) 
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CflLL  QRTDI(XST,YST,L1,L2) 

Q5LL  CFTDI(XFIN,YFIN,L3,L4) 

CPLL  a_INE(Ll,L2,Ll,L4,2) 

Q3U_  aLir«E(Ll,L4,L3,L4,2) 

CPLL  a_I^E(L5,L4,L3,L2,2) 

CflLL  QLirE(L3,L2,Ll,L2,2) 

C  TVC  FIRST  CLRV^  IS  DRPMN.  TTE  LlhE  CDLDR  IS  CHANGED  AM)  TVE 
C  SEECN)  OJFNE  IS  DRAkN 

CALL  (JTAa_(l,N5,XDAT,YDAT) 

IKOLCR  =  1 

CALL  aSETli=(MXJTS,IkrlJ3?,ISy^01.,KLRSV^1) 

CALL  OTABL(l,Nfe,XDATl,YDATl) 

CALL  OBEEP 
M  =  M  ♦  1 

C  IF  TTE  CLfIVES  APE  CF  DERIVATIVES  AND  A  LINEAR  FEBIDN  HAS  AJ  fEADY 
C  BEEN  SELECTED,  TFE  USER  MAY  DECIDE  TO  GE^ETV^TE  TVE  TAFEL  CDNSTWT5 
C  BASED  ON  THIS  PCRTIGN  CF  TIE  CLFVE. 

300  IF  ((  M  .GT,  1)  .AM).  (N4  .GE.O) )  TTEN 

CUES  =  GETERATE  TAFEL  CCNSTANT' 

CALL  □GTXT(2B,aLES,3,l,335,0) 

OLES  =  •  Y  OR  N  ?' 

CALL  aGrrXT(2B,aES,3,l, 320,0) 

CALL  CFEEP 
CALL  aiM<EY(I,l«B) 

IF  ((re  .EQ.  09)  .DR.  (MB  .ED.  121))  THEN 

flFS  =  '  • 

CALL  aGTXT(2B,aJES,3,l,335,0) 

CALL  GETXT(2B,aJEB,3,l,320,0) 

381  DUES  =  '  IS  THIS  TTE  AMDIC  DR' 

CALL  aSTXT(2B,aLES,3,l, 335,0) 

OLES  =  GyVTTCDIC  BRANCH' 

CALL  GETXT(2B,aLES,3,l,320,0) 
aES  =  '  A  DR  C  ?' 

CALL  aSTXT(2B,aLES,3,l,305,0) 

CALL  GBEEP 
C^EL  aiM<EY(I,l>B) 

IF  ((MB  .ED.  65)  .DR.  (MB  ,ED.97))  TVEN 
A1  =  KB 
A2  =  107 
GOTO  390 
ENDIF 

IF  ((TB  .ED.  67)  .CIR.  (TB  .ED.  99))  TTEN 
Cl  =  K6 
C2  =  K7 
GOTO  390 
ELSE 

com  3B1 
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EM)IF 

EMJIF 

IF  ((MB  .ED.  78)  .CR.  (MB  .ED.  IIO) )  THEN 
QDTO  390 
ELSE 

QDTO  380 
EM)IF 
bM>IF 

C  BY  TVE  USE  CF  TVE  (IFBCF  KEYS  ON  TVE  NUTERIC  KEYPW)  TVE  USER 

C  SELECTS  A  NEW  SET  CF  AXIS  TO  BE  FUJI  TED 

390  DIES  =  ■  ' 

CALL  (JGTXT(2B,C1JES,3, 1,335,0) 

CALL  aSTXT(2B,aLES,3, 1,320,0) 

C7)LL  aOTXT(28,QLES,3,l,305,0) 

GLES  =  KEW  AXIS  Y  OR  N  ?' 

CALL  CBTXT(2B,aL£S,3, 1,335,0) 

CALL  OHEEP 
CALL  QIN<EY(I,f>e) 

IF  ((MB  -ED.  B9)  .CF.  (MB  .ED.  121))  TVEN 
GOTO  400 
ELSE 

IF  ((MB  .ED,  78)  .OR.  (MB  .ED.  IIO))  TVEN 
GOTO  BOO 
FI  SF 

GOTO  390 
ENDIF 
EM)IF 

400  DIES  =  '  ■ 

CALL  aGrXT(2B,aLES,l, 1,335,0) 

CALL  (ETXT(2B,aUES,  1,1,320,0) 

GLES  =  'NEW,LAST,OR  ORIGINAL  (N,L,0) ' 

CALL  aGTXT(2B,QLES,3, 1,335,0) 

CALL  rw-y-p 
CALL  OIM<EY(I,MB) 

IF  ((re  .ED.  79)  .CR.  (MB  .ED.  Ill))  TTEN 
DO  410  I  =  1,N3 
XDAT(I)  =  TX(I) 

YDAT(I)  =  TY(I) 

410  aWTINLE 

DO  420  1  =  l,N4 
XDATl(I)  =  TXKI) 

YDATKI)  =  TYl(I) 

420  (DONTIIVLE 

M  =  O 
GOTO  5 
EM)IF 

IF  ((MB  .ED.  76)  .OR.  (MB  .ED.  lOB))  TFEN 
X5T  =  XSTl 
XFIN  =  XFINl 
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GOTD  530 
EM)  IF 

IF  ((  MB  .EQ.  7B)  .OR.  (MB  .EQ.  IIO) )  TVEN 
xsTi  =  xsr 
VSTl  =  YST 
XFINI  *  XFIN 
YFINl  =  YFIN 
M5  =  N5 
M6  =  N6 

OUEB  =  'USE  7FE  CXFGCF?  KEY  TO  MCVE' 

CflLl.  aBTXT(2B,(UES,3, 1,335,0) 

QLES  =  -70  TIC  LOWER  LEFT  HPM)' 

CflLL  aGTXT(2B,aLES,3,l,320,0) 

QUES  =  OORCR,  HLTE  TO  EXIT' 

OflLL  aGTXT(2B,QLES,3,l,305,0) 

CflLL  nW-FP 
CflLL  QLSTTATd) 

CALL  (l.TVPE(3,3) 

CflLL  ajtW(3aO,175) 

CflLL  Q1M<EY(I,L) 

IF  (L  .EQ.  72)  T>CN 
CflLL  GLMMT(0,5) 

GOTO  500 
ENDIF 

IF  (L  .ED.  75)  TWEN 
CALL  GLWMT(-5,0) 

GOTO  500 
EM)IF 

IF  (L  .EQ.  77)  TFEN 
CA_L  QLMDW(5,0) 

GDTO  500 
EMJIF 

IF  (L  ,ED.  80)  TVEN 
CA_L  GLr«MT(0,-5) 

QDTD  500 
EM)IF 

IF  (L  .ED.  71)  TVEN 
CA_L  GUCW(XL,YL) 

CAJ_  GLMCMT(320-XL,175-YL) 

CA-L  aiTOR(XL,YL,XST,VST) 

X2(l)  =  XST 
Y2(l)  =  YST 
CAJ_  DBETIjPI 5,3,43,3) 

CPU-  (JTABL(0,N7,X2,Y2) 

GOTO  510 
ELSE 
QUID  430 
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EM)IF 

GLES  =  UBE  7>t  OJRSDR  KEY  TO  l“0>E' 

CflLL  aGrrXT(28,GLJES, 3, 1,335,0) 

OLES  =  TO  T>E  LPPER  RIQTT  HAM)' 

CALL  aGTXT(2B,QL£S,3,l, 320,0) 

CUES  =  CDRVER,  HDTt  TO  EXIT- 
CALL  aETXT(2B,aLES, 3, 1,305,0) 

CALL  OBPPP 
CALL  QIM<EY(I,L) 

IF  (L  .ED.  72)  TVEN 
CALL  ajO)R(0,5) 

GOTO  520 
EM)IF 

IF  (L  .ED.  75)  TTEN 
CALL  a_m«(-5,0) 

GOTO  520 
EMJIF 

IF  (L  .ED.  77)  Tl-EN 
CALL  QLrtMR(5,0) 

GOTO  520 
EM)IF 

IF  (L  .ED.  80)  TTEN 
CALL  0LrCMR(O,-5) 

GOTO  520 
EM)IF 

IF  (L  .ED.  71)  TTEN 
CALL  ajTEAD(XR,YR) 

CALL  aLMGW(320-XR,175-VR) 

CALL  aiTOR(XR,YR,XFIN,YFIN) 

X2(l)  =  XFIN 
Y2(l)  =  YFIN 
CALL  aTABL(0,N7,X2,Y2) 

CALL  QLSTATIO) 

GOTO  530 
ELSE 
GOTO  510 
EM)IF 
EM)IF 

IF  (VNATE  .ED.  ' VOLTS' )  TVEN 
GOTO  780 
EM)IF 

K6  =  1 
K7  =  N3 

DO  600  I  =  1,N3 
X  =  TX(I) 

Y=TY(  I ) 

IF  (((X  .GT.  XST)  .AM).  (X  .LT.  XFIN))  .AM).  ((Y  .GT.  YSf) 
.AM).  (Y  .LT.  YFIN)))  TVEN 
K6  =  I 
GOTO  605 
ENDIF 


CONTINUE 
DO  610  I  =  K6,N3 
X  =  TX(I) 

Y  =  TY(I) 

IF  (((X  .LT.  XBT)  .OR.  (X  -GT,  XFIN) )  .OR.  ((Y  .LT.  YST) 
-OR.  (Y  .GT.  YFIN)))  TFEN 
K7  =  I 
GOTO  630 
EM)IF 
CDNTINLE 
N5  =  K7-K6 
DO  670  I  =  K:6,K7 
YDATI  I-K6*-l )  =  TY(I) 

XDAT(I-K64-1)  =  TX(I) 

CavfTIMJE 
KO  =  1 
K9  =  N4 
K2  =  N4 
IC3  =  1 

DO  700  I  =  1,N4 
X  =  TXl(I) 

Y  =  TYld) 

IF  (((X  .GT.  XST)  .«SO.  (X  .LT.  XFIN))  .AM).  ((Y  .GT.  YST) 
.AM).  (Y  .LT.  YFIN)))  TVEN 
KB  =  I 
GOTO  715 
EM)IF 
CCKTIMJE 
DO  720  I  =  KB,N4 
X  =  TXl(I) 

Y  =  TYl(I) 

IF  (((X  .LT.  XST)  .OR.  (X  .GT.  XFIN))  .OR.  ((Y  .LT.  YST) 
.OR.  (Y  .GT.  YFIN)))  TFEN 
K9  =  I 
GOTO  730 
EM)IF 

CDNTIMJE 
DO  770  I  =  KB,K9 
YDATld-KB*^!)  =  TYld) 

XDATld-KB^-l)  =  TXld) 

CONTINLE 
N6  =  K9  -  KB 
GOTO  225 
N5  =  O 
N6  =  O 

DO  785  I  =  1,N3 
X  =  TXd) 

Y=TYd) 

IF  (((X  .GT.  XST)  .AM).  (X  .LT.  XFIN))  .AM).  ( (Y  .GT.  YST) 


+  (Y  .LT.  YFIN)))  THEN 

N5  =  N5  +l 
XDAT(N5)  =  TX(I) 

YDAT(N5)  =  TY(I) 

EM)  IF 
CINTINLE 
DO  790  I  =  1,N4 
X  =  TXl(I) 

Y=TY1(I) 

IF  (((X  .GT.  XST)  .AM).  (X  .LT.  XFIN) )  .PND.  ((Y  .CT. 
+  .AM).  (Y  .LT.  YFIN)))  T)-EN 

N6  =  N6  +1 
XDATKM.)  =  TXKI) 

YDAT1(N6)  =  TYl(I) 

EM)IF 
CCmiNLE 
ri  =  O 
GOTO  225 

GLEB  =  DO  YOU  WANT  TO  RUT' 

CALL  QGTXT(2B,aJE5,3,l,335,0) 

GLES  =  THIS  Q¥EH  Y  OR  N  ?' 

CALL  GETXT(2B,ajES,3, 1,320,0) 

CALL  QEEEP 
CALL  QIM<EY(I,M8) 

IF  ((hB  .ED.  09)  .CF.  (re  .ED.  109))  TVEN 
GOTO  030 
ELSE 

IF  ((MB  .ED.  70)  .OR.  (MB  .EQ.  IIO) )  TVEN 
GOTO  050 
ELSE 

GOTO  000 
EM)IF 
EM)IF 

CALL  FUJT(XDAT,YDAT,N5,XDATl,YDATl,M>,XNArE,YNPrE, TITLE, 
+XMAX ,  XMIN,YMAX ,  YMIN,  XST ,  XFIN,  Y5T ,  YFIN,  XMftJCK.VMAJCR ,  TITLE2, 
-»FACT,NPrE) 

IF  (YNATE;  .ED.  VOLTS'  )  TVEN 
GOTO  390 
ENDIF 

IF  (YNATE  .EQ.  Dv/Di  )  TVEN 
GOTO  300 
ENDIF 

DO  860  I  =  1,20CO 
TX(I)=0 
TY( I )=0 
TXl(I)  =  O 
TYl(I)  =  O 
XDAT(I)  =  O 
XDATl ( I )  =  O 
YDAT(I)  =  O 
YDATKI)  =  O 
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THIS  SLERXJTIAE  R-OT5  T>£  9¥iFH5  AS  SHOkM  ON  TVE  SCREEN  DURING 
THE  CflLLINB  SLBRDUTIt^  GRPPHl. 

THE  MAJORITY  OF  TVE  PRDGRPM  IS  I^MX.VED  IN  SETTING  UP  flRGLJE^ENT5 
TD  BE  PASSED  TO  SUBRO/Tire  C□N^AI^ED  WITHIN  TIE  FUTTMATICS 
LIBRPRY.  THIS  MPNJAL  SHOULD  BE  CONSULTED  IN  ORDER  TD  DETERHIE 
TIE  USE  OF  TIE  flRBL£^ENTS  OTTER  THflN  TKBE  DEFI^ED  BELOW 

XDAT  :  TTE  ARRAY  C0N^AINI^e  TKE  X  VALUES  OF  THE  FIRST  CURVE 
TD  BE  RJJTTED. 

YDAT  :  TIE  FIRST  CIFVE'S  Y  S«?LUES 

XDATl  :  TIE  ARRAY  CCNATININB  THE  Y  VALUES  OF  THE  SECOND  CURVE 
TO  BE  RjOTTED. 

YDATl  :  TIE  SECOM)  CURVE'S  Y  VAUJES 
TITLE  :  TIE  TITLE  OF  TIE  PLOT 
TITLE2:  AN  DISORIBTION  OF  TFE  GRAFH 
XNPTE  :  TIE  X  AXIS  TITLE 

SATEUE:  TIE  Y  AXIS  TITLE  RCPERLY  POSITION  AM)  CATENTATED 

SUEROLITIIE  FUjr(XDAT,YDAT,M^l,XDATl,YDATl,hPT2,XNAIE,VNArE 
+, TITLE,  XMAX ,  XMIN,YMAX  ,VMIN,  XST  ,IEIN,  VST  ,YFIN,  XMAJOR,YMAJOR 
+,T1TLE2,FACT,  NA^E) 

REAL*4  XDAT(2000),YDAT(2CXX»,XDAT1(20CX»,YDAT1(2000) 

1NTEGER*2  JCDL(IO) ,JRDW(10) ,IET1,NFT2 
INTEGER  13, 14, 15,16,17, 18,111, 

CHAHACTER»8  NAfE,VNArE*5,XNAME*ll,TnL£*15,SArELE»12,TITLE2t32 
FACT*6,A4«1,  aLES«2B 
A4  *  '  ' 

SATELE  =  YNArE//A4//FAi:T 

XORG  =  XST 

VORG  =  YST 

DTX  =  AB5(XMAX-XMIN) 

DX  =  DTX  /MAX(^FT1,^PT2) 

DTY  =  ABS(YnAX-YMIN) 

lOPT  =  O 

r«<NT  =  5 

MXJTS  =  O 

LABEL  =  -1 

MEE  =  3 

ISIZE  =  1 

ISTEP  =  2 

IPORT  =  1 

CALL  ZINIT(  ISIZE,  ISTEP,  IPORT,BATCH) 

JCL  =  20 
JCR  =  620 
ICL  =  10 
ICR  =  320 

13  =  12 

14  =  150 


I 


CPU.  zaj*/(jcL,ia_,i^») 

CALL  ZCDr>W(JCR,ICR,16) 

CflLL  ZCaW(I3,I4,16) 

C«JL  ZCDW(  15, 16,16) 

CflU.  ZCIIVf/(  17,18,16) 

CRLL  ZCnvW(  Ill,  112,16) 

XMINl  =XmN  -.2*DTX 

xmxi  =  xmx+  .i$(dtx) 

YMAXl  =  vmx  +  .1*DTY 
YMINl  =  VlilN  -  .1»DTY 
rCEC  =  2 

C«_L  ZHjaT(JCL,JCR,ICL,ICR,XI1INl,XMftXl, YMINl, YWXl,XaRG,YtF«) 
C«_L  ZRTDI(XST,YST,L1,L2) 

CflLL  ZRTDI(XFIN,VFIN,L3,L4) 

C«_L  ZSETIP(-1, 1,0,1) 

MINOR  =  1 

CPUL  ZXAXIS(XSr,XFIN,)m]CR,MI^CR,L/«EL,^CEC) 

CflLL  ZYAXIS(YST,YFIN,YMftJ0R,MINCP,L/«EL,M3EC) 

C«_L  ZXTICS(10,JCa_,JRDW) 

CPLL  ZYTIGS(10,JCa_,JRCW) 

CflLL  ZLIf«(U.,L2,Ll,L4,l) 

CALL  ZL1ME(U.,L4,L3,L4,1) 
cm.  ZLirc:(L3,L4,L3,L2,l) 

cm  ZLI^E(L3,L2,L1,L2,1) 
cm.  ZPU 
OUES  =  '  ' 

cm  aErrxT(2B,ajEB,3,i,335,o) 
cm_  a8TXr(2B,QLES,3,l,320,0) 
cm.  a8TXT(2B,aLJE5,3,l,305,0) 

□LES  =  '  IF  EE5IFED  INGERT  A  ^EW  PEN' 

CALL  aGTXT(2B,aLES, 3, 1,335,0) 

CUES  =  'TFEN  FfESS  ANY  KEY' 
tm_  aGTXT(2B,aLJES,3,l, 320,0) 

CPLL  ntVFP 
CALL  QIM<EY(I,MB) 

CPLL  ZTPBL(l,NPri,XDAT,YDAT) 

CPLL  ZPU 

CUES  =  '  IF  DESIRED  INSERT  A  FEW  PEN' 

CPLL  Gl3TXT(28,aLES,3,l,335,0) 

CUES  =  'TIEN  PIESS  PNY  KEY' 
cm  aGTXT(2B,aLES,3, 1,320,0) 
cm  GBEEP 
CPLL  aiM<EY(I,r«) 

CPLL  ZTPBL(1,W^2,XDAT1,YDAT1) 

CPLL  ZPU 


cflLL  aiM<Ey(i,m) 

ni=ft  =  •  IF  DESIRED  IfGERT  A  ^EW  PEN' 
CPLL  QGTXT(2B,aLES,3,l,3:S»0) 

CUES  =  'TI-EN  PRESS  ANY  KEY' 

CALL  QBTXT(2B,aLES,3,l, 320,0) 

OUL 

CALL  ZPTXTA(ll,XNflrt,l) 

CALL  ZSETUP(2,1,0,1) 

CALL  ZIW(  1,1, 1760, 1400,0) 

CALL  ZPTXTI 15, TITLE,  1,15, 16) 

CALL  ZPTXT(8,NPfC,l,Ill,I6) 

CALL  ZPTXT(32,TI7LE2,1,I5,I8) 

CALL  ZEXTXT(0,0,0,0,0,1) 

CALL  ZPU 

CALL  ZPA(I3,I4) 

call  ZLP«EL(12,SflrPLE) 

CALL  ZFINIS 
END 
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THIS  IS  TVt  FIRST  SUBROUTINE  CPU  ED  U^EN  USING  THE  CUBIC 
SPLIINE  rETHDD. 


rDATA 

XDRT 


XDATl 

YDATl 

DIV 


TVC  PIRRAY  CCNTAININB  POTENTIAL  AM)  LCE  CLFFENT  r)E^EITY 
THE  AHTAY  CONTAINING  TVE  LOG  CURRENT  VALLES  FOR  TTE 
CATHODIC  BRANCH 

TFE  AFWAY  CCNTAININB  TVE  DERIVATIVES  FDR  TTE  CATHODIC 
BRANCH 

SAME  A  XDAT  BUT  FDR  TTE  ANODIC  BRANCH 

S»E  AS  YDAT  BUT  FCF  TFE  AMDIC  BRATCH 

TVE  ARRAY  CDNTAINIFG  TVE  DERIVATIVES  l«HICH  IS  BEING 

RETLIN  FRCn  SLBSEQLENT  CALLED  SUBROUTIVES 


TVE  INTEGER  V»ARI#»-ES  ACTUALLY  UEED  IN  THIS  SUBRDUTIVE  ARE 
USED  TO  DEFIVE  TVE  REGIOV6  WITHIN  BDTIH  BRAMVES  AVAILABLE 
FOR  USE 


SLBRDUTI^E  SLCPE1(MDATA,N,  3,NArE,Al,PO,Cl,C2) 
fEAL*4  XDAT ( 2000 ), YDAT ( 2000), XDATK 2000), YDATl (2000) 
FEAL»8  M)ATA(2000,2),DIV(1000) 

INTEGER*2  I  ,J  ,N,K,MI  ,N3,N4,N1,K1,A1,A2,C1,C2 
CHARACTER  NAfE*a,YNArE<5,TITLE2*32,TITLE3*50 
N1  =  N 

YNATE  =  Dv/Di' 

TIT1JE2  =  TAPEL  SLOPES:  CUBIC  SPLIIE' 

TITLE3=’  ' 

TVE  CATHODIC  BRAMH  IS  OEEKED  FOR  AN  INCREASE  IN  LOG 
CURFENT  DEVeiTY  WHICH  CDLLD  C^UEE  TVE  CALOJLATICN  OF  A 
DERIV#)TIVE  CF  IVFINTE  SLOPE 

DO  20  I  =  J-1,2,-1 

IF  (M)ATA( 1-1,2)  .LE.  MDATA(I,2))  GOTO  25 
ONTINUE 

K  =  riAX((J-250),l) 

GOTO  26 
K  =  I 
MI  =  J-9 
J1  =  K 

N3  =  ABS(K-MI) 
n  =  7 

CALL  CSR-IN(M)ATA,N,M,DIV,K,N3) 

DO  100  I  =  l,N5-3 
XDAT(I)  =  M)ATA(J-4-I,2) 

YDAT(I)  =  DIV(N5-I+1) 

CDNTIME 
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T>e  ANODIC  BRANCH  IS  CHECKED  FOR  LIXJ  CURRENT  VALUES  WHICH 
COULD  CAUSE  TEE  CALCULATION  OF  A  DERIVATIVE  CF  IIFINITE  SLIPE 


DO  120  I  =  J+1,N1 

IF  (M)ATA(I,2)  -LE.  M>ATA( 1-1,2) )  GOTO  125 
CDNTIMJE 

K  =  MIN((J+250),(N1)) 

GOTO  126 
K  =  I 
J1  =  J 
MI  =  K  -  9 
m  =  I»B(MI-^) 

CALL  CSPLIN(M3ATA,N,M,DIV,J,N4) 

DO  200  I  =  1,N4 
XDATl(I)  =  «)ATA(Jl+I+4,2) 

YDATl(I)  =  DIV(I) 

CONTINUE 

CALL  GRAFHl(XDAT,YDAT,N3,J,NArt,XDATl,YDATl,N4,YNArt,TnLE2, 
+A1  ,A2,Cl,C2,TI'n-£3) 

EM) 
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THIS  SUBROJrir^  SETS  UP  THE  MftTHIX  CF  CIEFFICIEIvfTS  TD 
BE  SOLVED  IN  ORCEH  TD  DETERMIlse  TVC  DERIVATIVE  USING  THE 
CUBIC  SPLIKE  rETHDD. 

M3ATA  :  TVE  ARRAY  CONTAINIISC  POTENTIAL  LOG  CLRRENT  DENSITY 
CMAT  :  TEE  7x7  MATRIX  CONTAINING  THE  COEFFICIENTS 
YMAT  :  TFE  7x1  VECTOR  CONTAINING  TVE  RIGHT  HAM)  SIDE  OF 
TIE  SYSTEM  CF  EOJATIONB.  WEN  FETURED  FROM  TVE 
9C1.VING  SUBROLm^E  IT  CTFOAINB  TVE  SOLLfTICN. 


SLERC)Lm^E  CSPLI^E(^DATA,J,M,DIV,K,MI) 

REAL*8  XDAT(IOOO),  YDAT(IOOO) ,M)ATA(2CXX),2)  ,YMAT(7) 
fEAL«a  CMAT(7,7),AI,BI,CI 
fEAL<a  DlVdOOO) 

iNiEOEnr:’  n,mi,k,j 

DU  10  I  =  1,7 
DO  5  J  =  1,7 

CJ1AT(I,J)  =  O.DO 
5  CENTIME 

10  CONTINLE 
M  =  7 
K  =  K-1 


DO  200  N  =  1,MI 


CMAT(1,1)  =  2.D0*(MlAIA(K>2,2)-M)ArA(K,2)) 

CmT<l,2)  =  M)ATA(K+2,2)  -  MDATA(K+1,2) 

CriAT(7,7)  =  2.D0*(M3ATA<K+e,2)-MDATA(K4ilb,2)) 

CJ1AT(7,6)  =  M)ATA(K+e,2)-M)ATA(K+7,2) 

YMAT(1)=  ((MDATA(K+2,l)-MDATA(K+l,l))/(M3ATA(K+2,2)-M}ArA 
+  (K+1,2)))  -  ( (r«ATA(K+l,l)-ff)ATA(K,l))/(M)ATA(K+l,2)-M)ATA 
>(K,2))) 

YMAT(l)  =  YliAT(l)*6. 

YMAT(7)=  ((M)ATA(K+8,l)-M)ATA(K+7,l))/(rt>ATA(K40,2)-M)ATA 
f(K+7,2)))  -  ((M)ATA(K+7,l)-MDATA(K+<>,l))/(rDATA(K+7,2)-M)ATA 
>(K+6,2))) 

YMAT(7)  =  YMAT(7)*6, 

DO  100  I  =  2,6 
K  =  K  +1 

CMAT(I,I)  =  2.DO«(  (M)ATA(K+2,2)-M)ATA(K,2) ) ) 

CriAT(  1,1+1)  =  M)ATA(K+2,2)-M)ATA(K+l,2) 

CMAT( 1,1-1)  =  M)ATA(K+1,2)  -  M)ATA(K,2) 

YMAT(I)  =  ((M)ATA(K+2,1)-MXITA(K+1,1))/CMAT(  1,1+1)) 

+-( (VDATACK+l  ,1  )-Vt)ATA(K,l )  )/CMAr(  I ,  I-l ) ) 

VMAT(I)  =  YMAT(I)*6. 
lOO  CUfTlNJE 
K  =  K  -  5 

CALL  LIhEYl (CHAT, YMAT, M) 

AI  =  VMAT(5)  -  YMAT(4) 


ifdkllMilklii 
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AI  =  Ar/((6.D0)»((rt3ATA(KM,2)-rt3ATA(K+3,2)))) 

BI  =  YMAT(4)/2.D0 

Cl  =  (rt)ATA(K+4,l)-MDATA(K+3,l)  )/(l«CATA(K+4,2)-M)ATA(K+3,2) ) 
Cl  =  CI-<l/6)»(((2«VMAT(4))+YMAT(5))$(t«ATA(K+4,2)-t®ATA(K+^ 
+3,2)) ) 

DIV(N)  =  Cl 

CCNTINUE 

EM) 


THIS  SUBRCUTire  SU_'-ES  TfC  7x7  SYSTEM  CF  EOUftTlQNB  GE>EF!ATED 
BY  CSPLIN.FQR.  TVE  PROmPM  WPS  ORIGINPLLY  FFCM  THE  TEXTBDCK 
"FCFnrVlN  77  for  Enqineers  aid  Scientists”  BY  NYHDFF  PM) 

LEEHTMA,  PAGES  264-265.  TFE  MAIN  MDDIFICATICN  IS  TT-E  ELIMINATIC3M 
□F  TFE  VfiRI#«_E  DIFENBIGNING  CF  ATFWYS  RESENT  IN  TFE  CFIGINPL 
PRCCFWM. 


SGBR□L^■I^E  LINBV1(C3>1AT,YMAT,N) 

INTECEP  N, PIVOT 

FEflL»8  «JG(7,a),  CIMAT(7,7),  YMAT(7)  ,X(7)  ,MJLT,TOF» 
DO  20  I  =  1,  N 
DO  10  J  =  1,N 

«JG(I,J)  =  ODMATd.J) 

cumtime; 

COFTTINE 
DO  30  I  =  1,N 

AUB(I,N*^1)  =  YMAT(I) 

comtime; 

DO  70  I  =  1,N-1 

IF  (AIJG(I,I)  -EQ.  O.)  TFEN 
PIVOT  =  O 
J  =  I+l 

IF  ((PIVOT  .EQ.  O)  .PfC.  (J  .LE.  N) )  TFEN 
IF  (AUG(J,I)  .NE.  O.)  PIVOT  =  J 
J  =  J+1 
GOTO  35 
EM)IF 

IF  (PIVOT  .ED.  O)  TFEN 
RONT*,  SIFOULPR- 
STCP 
ELSE 

DO  40  J  =  l.Nd 
TETP  =  ALEd.J) 

PUGd.J)  =  ALE(PIVOT,J) 

«JG(PIVOT,J)  =  TEFF 
dNTINLE 
ENDIF 
EM)IF 

DO  60  J  =  I+1,N 

MET  =  -«JB(J,I)/«-Dd,I) 

DO  50  K  =  I,N«-1 

ALE(J,K)  =  «JB(J,K)  +  MJLT*ALE(I,K) 

CXFfTIME; 

ONTIIyLE 

CXNTIMJE 

X(N)  =  «JG(N,N*1)/PLJG(N,N) 

DO  90  J  =  N-1,1,-1 
X(J)  =  «JG(J,Nd) 
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m  00  K  =  JH,N 

X(J)  =  X(J)-(«JG(J,K)*X(K)) 
0O  tXlsfriMJE 

X(J)  =  X(.I)/eLG(J,J) 
so  CCNriNUE 

DO  100  I  =  1,N 
YMATd)  =  X(I) 
lOO  CONTir^ 

END 
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*  THIS  SLJERXfTI^E  CUvKyERTS  A  NUhERIC  VARIBLE  (X6)  CF  F0R1AT 

*  +/-  XX.XXXX  TQ  A  CHPRACTER  WKEE  WIDTH  IS  8  CHPRFCTERS. 

«  THIS  IS  hECEBSmy  AS  TFE  H-OTMATICS  LIBRflRY  IS  NOT  EQUIPPED  WITH 
«  A  SUBROJTINE  FCR  PLOTriNB  NUTtRIC  VARIBLEB. 


SUBROJriNE  CnsW(A7,X6) 

I^fIEGER  01,02,03,04,05,06 

CHPHACTERIl  A6(0:9) ,  C1,C2,C3,C4,C5,G:6,C:7,C8,A7*8 
A6(0)  =  O' 

A6<1)  =  I- 
A6(2)  =  '2' 

A6(3)  =  ■3- 
A6(4)  =  4' 

A6(5)  =  '5- 
A6(6)  =  6' 

A6(7)  =  ■7- 
<Vi(8)  =  B' 

#«.(9)  =  '9' 

Cl  =  '  • 

C2  =  '  ' 

C3  =  '  ' 

04=*' 

C5  =  •  ' 

C6  =  •  ’ 

C7  =  '  ' 

IF  (X6  .LT.  0.0)  Cl  = 

X6  =  #»S(X6) 

XI  =  X6/10. 

01  =  IFIX(Xl) 

X2  =  X^r-<10.*01) 

02  =  IFIX(X2) 

X3  =  X6  -  (10.401)  -  02 
X3  =  X3/.1 
03  =  IFIX(X3) 

X4  =  (X3-a3)/.l  +  .001 

04  =  IFIX(X4) 

X5  =  (X4-04)/.l 
05  =  IFIX(X5) 

X7  =  (X5-05)/.l 
06  =  IFIX(X7) 

DO  lO  I  =  0,9 

IF  (01  .ED.  I)  C2  =  A6(I) 

IF  (02  .ED.  1)  C3  =  ^(I) 

IF  (03  .ED.  I)  04  =  A6(I) 

IF  (04  .ED.  I)  CS  =  (¥j(I) 

IF  ((35  .ED.  I)  C6  =  A6(l) 
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THIS  SLBRDUrir>t  SETS  SSTE  OF  I>t  CHP«fiCTER  STRIM3S  LEED 
l»FEN  USING  7T-E  CENTRAL  DIFFEIEME:  rETHDD 

SUBHOJTirE  3LaFE(rf)ATA,N,J,NPrE,Al,A2,Cl,C2,Nl,N2,N3,N4) 
REPL*4  XDA T ( 2000 ), YCAT ( 2000 ), XDATK 2000 ),YDAT1( 2000) 

FEPL*8  MDATA( 2000,2) 

INTEGER*2  I  ,J  ,N,A1  ,A2,C1  ,C2 

CHWfiCTER  N«ME*8,YNArt*5,TITLE2*32,  LISrr*5O,TI7LE3*50 

YNRTt  =  Dv/Di'  _ 

riTLE2  =  TAFEL  SLOPES:  CENim.  DIFFEHiNCE' 

TITLE3=‘ 

Cm_  DATADEL(NDATA,J,YDAT,YnATl,XDAT,XDATl,N2,Nl) 

CALL  GRAPHl(XDAT,YDAT,N3,vI,NAft,XDATl,YDAri,N4,YNArt,TITL£2, 
+A 1 ,  A2 ,  Cl ,  C2 ,  T I  TU:3 ) 

EM) 


riwP»,i 


THIS  SUERUTire  COvOJCTS  MJr^RICAL  DIFFERENTIATION  LEINS 
TIE  4  POINT  CENTFJPL  DIFFEFNCE  ftTHDD. 

rOATA  =  PHRAY  CONTAINIF^  POTENTIAL  AM)  LOG  I  DENSITY 
SLCPE  =  AfffWY  CONTAINING  DE/D(LDGI )  AND  LOG  I 
FSLCFE  =  CENTFJAL  DIFFERCE  DIFFERENTIATION 
ECDFRl  =  LPPERLIMIT  CF  ANDDIC/CATHDIC  DISTRIBUTION 
ECCRR2  =  LCWER-IMIT  CF  ANODIC/CATFCDIC  DISTRIBUTION 
N  =  MJTOER  OF  DATA  POINTS 

J  =  LOCATION  WITHIN  FOATA  CF  ECOFR 

I  =  CATHODIC  BRPNCH  LOOP  COLNTER 

K  =  ANODIC  BRANCH  LOOP  COLNTER 

H  =  DIFFERENTIAL  VALLE  IN  POTENTIAL 

II  =  VALLE  CF  LOG  I  AT  E-2H 

12  =  VALLE  OF  LOG  I  AT  E-*-' 

13  =  VALLE  CF  LOG  I  AT  E 

14  =  VALLE  CF  LCB  I  AT  EHT 

15  =  VALLE  CF  LOG  I  AT  E+2H 

DEDI  =  SLOPE  VALLE  OF  13 

SLJBFaJTIfE  DATADEL  (M)ATA,J,YDAT,YDAT1,XDAT,XDAT1,N2,N1) 
FEAL»8  M)ATA<2CXX),2),X1,X2,X3,X4,X5,H,FSLDFE,DEDI 
FEAL«4,  XnATl(lCX)0),YDATl(lCXX)),  YDAT(  1(XX»  ,XDAr ( 1(X)0) 
INTEGERTZ  J,K,I 

FSLCFE  (H,Xl,X2,X4,X5)=(Xl+(a.DO*X4)-(8.DOtX2)-X5)/(12.DO*H) 
FEFFER1ING  CALCULCATICfe  ON  CATHODIC  BRANCH 
K  =  O 

DO  100  I  =  J-2,Nl+2,-l 
K  =  K>1 

H  =  MDATA(I,1)  -  MDATA(I-1,1) 

XI  =  rt)ATA( 1-2,2) 

X2  =  MDATA(  1-1,2) 

X3  =  M)ATA( 1,2) 

X4  =  M)ATA( 1+1,2) 

X5  =  rCATAC 1+2,2) 

DEDI  =  l.DO/FSLCFE(H,Xl,X2,X4,X5) 

YDAT(K)  =  DEDI 
XDAT(K)  =  ^CATA(I,2) 

CONTINLE 
I  =  O 

PERFDFMING  CALCLLATIONB  ON  TVE  AMBIC  ffWCH 

DO  300  K  =  J+2,  N2-2 
1  =  1  +  1 

H  =  M)ATA(K+1,1)  -  M)ATA(K,1) 


Xi  =  rt»TA(K-2,2) 

X2  =  ri»TA(K-l,?) 

X3  =  rt>ATA(K,2) 

X4  =  rt)ATA(K^l,2) 

X5  =  rt3ATA(K+2,2) 

CEDI  =  l.DO/FBLXFE(H,Xl,X2,X4,X5) 
YDATl(I)  =  CEDI 
XDATl(I)  =  r«ATA(K,2) 

CONTINUE 
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